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Abstract 


A need exists to safely integrate Unmanned Aircraft Systems (UAS) into the National 
Airspace System. Replacing manned aircraft’s see-and-avoid capability in the absence of an 
onboard pilot is one of the key challenges associated with safe integration. Sense-and-avoid 
(SAA) systems will have to achieve yet-to-be-determined required separation distances for a 
wide range of encounters. They will also need to account for the maneuver performance of the 
UAS they are paired with. The work described in this paper is aimed at developing an 
understanding of the trade space between UAS maneuver performance and SAA system 
performance requirements. An assessment of current manned and unmanned aircraft 
performance was used to establish potential UAS performance test matrix bounds. Then, near- 
term UAS integration work was used to narrow down the scope. A simulator was developed with 
sufficient fidelity to assess SAA system performance requirements for a wide range of 
encounters. The simulator generates closest-point-of-approach (CPA) data from the wide range 
of UAS performance models maneuvering against a single intruder with various encounter 
geometries. The simulator is described herein and has both a graphical user interface and batch 
interface to support detailed analysis of individual UAS encounters and macro analysis of a very 
large set of UAS and encounter models, respectively. Results from the simulator using 
approximate performance data from a well-known manned aircraft is presented to provide 
insight into the problem and as verification and validation of the simulator. Analysis of climb, 
descent, and level turn maneuvers to avoid a collision is presented. Noting the diversity of 
backgrounds in the UAS community, a description of the UAS aerodynamic and propulsive 
design and performance parameters is included. Initial attempts to model the results made it 
clear that developing maneuver performance groups is required. Discussion of the performance 
groups developed and how to know in which group an aircraft belongs for a given flight 
condition and encounter is included. The groups are specific to airplane, flight condition, and 
encounter, rather than airplane-only specific. Results and methodology for developing UAS 
maneuver performance requirements are presented for each maneuver as well. Results for the 
vertical maneuver indicate that a minimum specific excess power value can assure a minimum 
CPA for a given time-to-go prediction. However, smaller values of specific excess power may 
achieve or exceed the same CPA if the UAS has sufficient speed to trade for altitude. Level turn 
results are less impacted by specific excess power and are presented as a function of turn rate. 
The effect of altitude is also discussed for the turns. Next steps and future work are discussed. 
Future studies will lead to better quantification of the preliminary results and cover the 
remainder of the proposed test matrix. It is anticipated that this will be done in conjunction with 
RTCA SC-228 over the next few months. 
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1 Introduction 


The United States National Airspace System (NAS) has evolved to both support and rely on 
the capabilities of manned aircraft. Currently, the Department of Defense (DOD), other public 
entities, and private industry are showing increasing interest in gaining regular access to fly 
unmanned aircraft systems (UAS) in the NAS. One of the major barriers to integrating UAS in 
the NAS is the requirement to see and avoid other aircraft, a task currently performed by the 
onboard pilot. Title 14 of the Code of Federal Regulations (14 CFR) Part 91.113 [1] states, 
“vigilance shall be maintained by each person operating an aircraft so as to see and avoid other 
aircraft.” Thus, pilots are expected to see and avoid other aircraft and to maneuver in 
predictable ways to preserve the safety, orderliness, and efficiency of the NAS environment. 

UAS will be expected to operate in a similar manner, but with sense and avoid (SAA) systems 
replacing the see-and-avoid capability of manned aircraft. However, the acceptable design 
space and required capabilities for SAA in this environment remains largely undefined. 

According to the Second Caucus Report from the Federal Aviation Administration (FAA) 
sponsored Sense and Avoid workshop, “for a technical system to perform the function of a pilot 
to remain well clear, it is necessary to have an unambiguous, implementable definition of the 
separation minima” [2]. Without a quantitative definition of “well clear,” it is difficult to develop a 
UAS that is able to conform to current airspace rules. In addition to quantification of well-clear 
separation minima, understanding the trade space between UAS maneuver performance 
capabilities and SAA system requirements necessary to maintain the minimum separation 
standards is needed. Some separation boundaries are currently defined as the last time to 
maneuver to avoid penetration of the next boundary. This type of definition has UAS 
performance inherent in the definition and incorporates encounter parameters as well. Under 
these definitions, UAS performance and the well-clear or near mid-air collision (NMAC) 
boundaries must be combined to determine when the SAA system must initiate an avoidance 
maneuver. Therefore, as the FAA works to define new rules that will lay out the requirements for 
all UAS SAA systems, it is essential to understand and quantify the trade space between aircraft 
maneuver performance and the capabilities of the SAA system. UAS designers will also benefit 
from a deeper understanding of this trade space. 

A large body of work related to collision avoidance was generated during the development 
of the Traffic Alert and Collision Avoidance System (TCAS) [3]. TCAS is now required 
equipment on large manned aircraft acting as a last line of defense against collisions with other 
aircraft. TCAS analyzes the projected flight path of ownship and intruders and issues advisories 
to pilots based on a complex algorithm, thereby providing a layer of protection against collisions 
other than air traffic control (ATC) and the pilot’s vision and discernment. However, TCAS was 
designed and intended for use on aircraft having a large number of passenger seats and 
meeting standards defined in CFR 14 Part 25. This generally defines a class of aircraft that flies 
at high altitudes and at relatively fast speeds (transport category airplanes) when compared to 
many general aviation (GA) aircraft and UAS. TCAS relies on timely pilot responses and aircraft 
performance parameters to achieve collision-avoidance minimums based on “Resolution 
Advisory” times prior to a predicted collision. Transponders provide limited relative state 
information of intruder aircraft used by the threat detection and maneuver advisory algorithms. 
Technological limitations of transponders at the time TCAS was developed further restricted 
TCAS design options. Additionally, TCAS is intended for collision avoidance only and, therefore, 
does not address other separation requirements, such as remaining well clear as described 
earlier. As a result of these limitations, TCAS requirements do not take advantage of currently 
available sensor capabilities and require higher climb and descent performance than the broad 
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array of current manned aircraft exhibit. The widely diverse UAS fleet often exhibits even less 
maneuver performance currently, sometimes because of airframe and power-plant limitations 
and sometimes because of the flight-control system. 

These observations led to the development of a simplified simulation tool and design of a 
simulation experiment to analyze the interaction of UAS and SAA system performance 
requirements. The simulation tool allows broad and rapid variation of UAS design parameters 
that impact maneuver performance. For each simulation run, a UAS and a non-maneuvering 
intruder are set on a collision course, and the UAS immediately initiates a maneuver to avoid 
the intruder, with resulting Closest Point of Approach (CPA) and other airplane states recorded. 
The avoidance maneuvers evaluated were commanded climbs, descents, and level turns 
across a range of rates with airplane performance limitations impacting the resulting maneuver. 
Intruder azimuths in the UAS’ forward hemisphere were evaluated as well as climbing and 
descending intruders. The encounter scenarios resulted in CPAs spanning from NMAC ranges 
to well beyond those that could be considered well clear. Three simulation models were 
developed to account for UAS maneuver performance limitations associated with climb, 
descent, and level turn maneuvers independently. This approach significantly reduces model 
complexity and, therefore, the amount of time required for each simulation run. 

A wide range of diverse existing fixed-wing airplane types were used to define parameter 
ranges representing the broad spectrum of UAS with present and future aerodynamic and 
propulsive performance and maneuverability. The UAS parameter ranges were then narrowed 
to focus the simulation on performance ranges that could potentially limit conduct of avoidance 
maneuvers and to align with near-term UAS integration concepts as defined in the RTCA 
Special Committee - Minimum Operational Performance Standards for Unmanned Aircraft 
Systems (SC-228) Terms of Reference [4] and work plan [5], Even with the parameter range 
reductions, millions of combinations of airplane performance and encounter parameters were 
analyzed in an effort to describe the trade space between SAA system and UAS maneuver 
performance requirements. 

Design parameters that have first-order effects on airplane performance were used to 
develop the simulation test matrix resulting in a broad set of airplane maneuver performance 
capabilities. The parameter matrix was also used to determine likely initial UAS speeds. Rather 
than simulating specific airplanes, the approach was used to cover a broad range of maneuver 
performance. Before data gathering began, a verification and validation study was performed to 
ensure that the simulation accurately modeled maneuver performance across the parameter 
ranges of interest. Example results from one manned aircraft that was modeled are included. 

Using results from the simulation, a number of candidate performance parameters were 
evaluated to determine the best universal predictors of CPA. The objective was to identify a set 
of performance parameters that could be used to predict the avoidance maneuver result given 
current or future UAS. The performance parameters could also be used to identify when a 
maneuver must be initiated to generate a desired CPA from a series of encounters. To be 
useful, the performance parameters must be known or easily estimated for a given UAS or 
future UAS design. Extensive analysis was done of the vertical climb maneuver and a proposed 
approach for evaluating performance parameters was developed for aircraft that can sustain a 
commanded climb and those that cannot. A more limited analysis of the vertical descent and 
level turn maneuvers was conducted with lessons learned and methodology developed included 
in the paper. The descent maneuver is evaluated in detail in a companion document [6]. 
Airplanes that can sustain the commanded maneuver were analyzed for all three maneuvers. 
Approaches for analyzing the aircraft that cannot sustain a commanded maneuver are also 
discussed. Conclusions from analysis of performance-limited groups are for more specific sets 
of parameters with the focus on methodology. 

This paper shows results and trends in CPA associated with the best performance 
predictors. As indicated by the able to sustain, and unable to sustain, commanded maneuvers 
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references above, maneuver performance groups that are analyzed separately to make the 
problem more manageable became necessary and are discussed. The performance groups are 
specific to the UAS airframe/power-plant, commanded maneuver, and the encounter. A means 
to identify into which group an airplane, flight condition, and encounter falls is defined. Analysis 
of results shows how specific performance parameters can be used to identify when a 
maneuver must be initiated to generate a certain CPA from a series of encounters. In some 
cases the results are broad and in other cases analyzed against a particular CPA. In the latter 
case, the existing data and approach described can be analyzed against other CPAs in the 
likely event that values defined by regulators differ from those used in the examples herein. 

The results and methodology presented in this paper will help regulators determine 
minimum maneuver performance requirements for UAS operating in the NAS with a SAA 
system. Results will also help determine the point in the SAA system timeline at which the 
maneuver must be initiated, whether the desired maneuver is maximum performance or 
something less. The resulting trade space is expected to support development of flexible SAA 
system design requirements, thus accommodating as many UAS types in the NAS as possible. 

2 Performance Model and Simulation 

This section describes the simulation in detail and includes: the range of UAS modeled; the 
simplified simulation tool developed to enable rapid evaluation of millions of combinations of 
aircraft and encounter parameters; and the commanded avoidance maneuvers modeled. A 
description of validation and verification of the simulation conducted to ensure the simulation 
accurately captured flight performance is included. The section starts with a derivation of 
parameter ranges to be used and description of the modeling technique to provide the reader an 
understanding of the simulation fidelity. 

2.1 Derivation of UAS Parameter Ranges 

Since the goal of the simulation was to evaluate UAS performance in conducting avoidance 
maneuvers across a wide range of UAS, the UAS maneuver performance was modeled by 
varying key airplane design parameters across applicable ranges. A survey of existing airplanes 
was done to identify applicable parameter ranges of UAS and manned fixed-wing airplanes. 
Most current UAS fly slowly and have limited maneuverability, while others have performance 
approaching that of airliners and corporate jets. UAS with higher performance and 
maneuverability are expected in the future, but these aircraft do not represent a limiting case; 
thus, along with fighter aircraft, they are absent from the sample. 

A broad sample of existing manned and unmanned airplanes is shown in Table 1 to 
represent the spanning set of airplane performance of interest. Key design parameters for these 
airplanes, which coarsely define the performance of each, were obtained from publically 
available sources. Some of these airplanes have cruise speeds on the order of 60 knots true 
airspeed (KTAS) while others cruise at 500 KTAS or more. Some can barely climb at 500 ft/min 
and others can sustain thousands of ft/min. Some in the table have stall speeds that are higher 
than the velocity not to exceed ( V NE ) of others in the sampling. Finally, some of these airplanes 
cannot reach 18,000 ft Mean Sea Level (MSL) (bottom of Class A airspace) and others cruise 
above 30,000 ft MSL. 
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Table 1 : Selected aircraft and publicly available design parameters. Some values were 

estimated from publically available data. 


# 

Aircraft 

T/W 

min 

T/W 

max 

W/S 

(lb/ft 2 ) 

min 

W/S 

(lb/ft 2 ) 

max 

AR 

e 

(-) 

\DJ max 

C|_max 

Aircraft type 

1 

Aerosonde 

0.07 

0.10 

4 

6 

10 

0.90 

20 

1.85 

Small UAS 

2 

Boeing 757-200 

0.31 

0.63 

65 

131 

8.0 

0.85 

18 

2.75 

Transport 

3 

C-17 

0.28 

0.57 

74 

154 

7.2 

? 

? 

3.54 

Cargo Transport 

4 

Cessna 172 

0.22 

0.38 

9.7 

15 

7.5 

0.80 

9 

2.38 

General Aviation 

5 

Cessna R182 

0.25 

0.35 

13 

18 

7.4 

0.80 

11 

2.27 

General Aviation 

6 

Embraer 

EMB-145 

0.16 

0.33 

46 

77 

7.8 

0.80 

15 

2.35 

Regional Jet 

7 

Fokker 
Dr. 1 

0.30 

0.44 

4.4 

6.4 

4.0 

0.65 

8 

1.24 

Tri plane 

8 

Global Hawk 

0.26 

0.77 

16 

48 

24 

0.98 

30 

1.21 

Large UAS 

9 

Ikhana 

0.19 

0.40 

19 

41 

17 

0.90 

30 

1.88 

Large UAS 

10 

Nimbus 4M 

0.10 

0.13 

6.2 

8.6 

39 

0.98 

60 

? 

Motor Glider 

11 

Piper 

PA-28- 180 

0.15 

0.26 

7.7 

15 

5.6 

0.76 

10 

1.99 

General Aviation 

12 

Shadow B 

0.32 

0.35 

12 

13 

11 

0.90 

13 

1.60 

Small UAS 

13 

Stemme S10 

0.10 

0.13 

7.0 

9.3 

28 

0.98 

50 

1.54 

Motor Glider 


Using the broad sample of airplanes shown in Table 1 as a baseline, and applying related 
subject matter expertise, a test range for each UAS design parameter was selected for the 
current study, as shown by the blue lines in Figure 1 . Note that these test ranges fall inside the 
parameter range for the airplanes in Table 1 . In order to reduce the number of combinatorial 
runs to a more manageable size, a decision was made to focus the simulation on UAS 
performance characteristics that would cause more difficulty in conducting effective avoidance 
maneuvers and more closely align with near-term needs. A key reason for doing this study is to 
assess airplanes that do not have the performance capability of those for which TCAS was 
originally designed. A high-performance UAS can easily conduct an avoidance maneuver, while 
a UAS that can marginally maintain its desired straight, level, unaccelerated flight condition is a 
limiting case of interest. Note that marginally able to maintain straight and level unaccelerated 
flight is a fairly normal flight condition for GA airplanes and for some types of UAS. This 
condition can be mission related; for example, a slow flight speed can be used for maximum 
endurance. “Marginally able” can also mean straight and level at maximum throttle and speed 
with no excess energy, meaning any climb or turn maneuver would result in a loss of speed. 
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Figure 1 : Test matrix of aircraft design parameters. Red dots are the values from the 

sample aircraft in Table 1 with the number corresponding to the aircraft number in 
the table. The blue bars show where the test matrix boundaries were set. 


The reduced matrix allowed results to be published earlier and they cover UAS performance in 
the near term as well as the focus area associated with RTCA SC-228 Terms of Reference [4]. 

Each design parameter is briefly described in the subsection below and rationale for the 
selected range is discussed. The simulation input is made in terms of these key design 
parameters, and they are related to key performance parameters within the simulator. Both the 
design and performance parameters used are readily available to airplane designers. Thus, the 
results can be readily applied to existing and new UAS designs. The parameters chosen are 
essentially for a “unit” airplane, which is to say that given the same set of parameters, the 
performance of the resulting airplane can be estimated without knowing the size of the airplane, 
i.e., a wing span of 6 ft or a wing span of 200 ft with the same set of parameters used here 
would have the same performance. This allows for comparison of airplane performance 
regardless of the physical size of the vehicle. Thus, an airplane designed to carry a 10-lb 
camera or a 100-lb camera at the same speed and altitude ranges could have very similar sets 
of design parameters because they could be essentially scale models of each other. In reality it 
is not quite this simple, as specific structural mechanics, aerodynamics, and engine 
performance may differ, but it is illustrative of the value of using these parameters. 
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2.2 Parameter Values for Simulated UAS 


This section gives a brief description of each of the UAS design parameters used in the 
study and describes the test range used in the simulation and the rationale for choosing that 
range. 

2.2.1 Thrust-to-Weight Ratio (^) 

Range: 0.1 to 0.4 (non-dimensional) 

The ratio of thrust to weight is a measure of the propulsion system’s thrust per unit weight of 
the airplane. The range of 0.1 to 0.4 was chosen because the vast majority of existing airplanes 

have ^ < 0.4 for at least part of their flight regime and, as can be seen from Table 1 most do not 
exceed it. Low ^ values are experienced when most airplanes reach their service ceiling and 

are near their maximum gross weight. UAS with higher ^ values are not the focus of this study 

as they will be able to outperform the airplanes of interest and do not represent a limiting case. 
Many of the configurations tested had sufficient thrust for the commanded climb maneuvers as 
will be seen in the results section. 

2.2.2 Wing Loading (^) 

Range: 1 to 25 (lb/ft 2 ) 

Wing loading is the airplane’s weight divided by the wing area, measured in pounds per 
square foot. A — range of 1 to 25 lb/ft 2 was chosen as the range of interest for this work 

because it is generally associated with single-engine CFR 14 Part 23 airplanes . It does not 
nominally cover airplanes that fly fast or at very high altitudes which are left for later studies. 

2.2.3 Aspect Ratio (AR) 

Range: 5 to 40 (non-dimensional) 

Aspect ratio is a fineness ratio of the wing computed as wing span divided by mean 
aerodynamic cord. High AR airplanes are characterized as having lower induced drag resulting 
in higher lift-to-drag ratios and higher aerodynamic efficiency. Conversely, low AR airplanes 
generally are capable of greater roll accelerations and rates but have lower aerodynamic 
efficiency. The range of interest for this study spans a wide range from unusually low aspect 
ratio airplanes to high aspect ratio motor gliders. It does not include the extreme AR values of 
some high-performance manned gliders. Note that the highest aspect ratios (AR) and maximum 

lift-to-drag ratios ((-) ) in Table 1 are from high-performance sailplanes. 

2.2.4 Maximum Lift/Drag (- ) 

D max 

Range: 5 to 50 (non-dimensional) 

1 Part 23 single engine land aircraft will generally have j less than 25 lb/ft 2 because their stall speed 
cannot exceed 61 KCAS per 14 CFR Part 23 [7] 
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This is the maximum lift-to-drag ratio for the airplane configuration. This value is used in the 
simulation to define the parasite drag of the airplane, a value that is difficult to obtain otherwise. 
The selected range of interest includes all airplanes in Table 1 but excludes a limited set of 
very-high-performance sailplanes. 

2.2.5 Maximum Coefficient of Lift 1C, ) 

Range: 0.6 to 2 (non-dimensional) 


The maximum coefficient of lift coupled with — is used to define the minimum velocity (stall 

speed) of the airplane configuration. The 2-D airfoil theory (potential flow) used in the simulation 
does not predict stall. Thus, C Lmax is an independent variable in the study, effectively imposing a 
stall constraint. It is also used in conjunction with initial percent of C Lmax to specify the airplane’s 
initial speed. This study evaluated airplanes at initial speeds ranging from nominal cruise to 
loiter. 

2.2.6 Initial Percent of Maximum Coefficient of Lift ( %C L ) 

Range: 0.08 to 0.80 (non-dimensional) 

This variable is used in conjunction with C, and — to define the initial velocity of the 

airplane. The approach sets up conditions associated with different phases of flight and parts of 
the flight envelope. It avoids direct association with angle of attack, which has different zero 
reference lines associated with specific configurations. As mentioned in the previous section, 
the range chosen spans from a nominal cruise condition to a loiter condition. Loiter is 
considered to be a best endurance condition, typical of many current UAS missions, and is 

therefore close to the - condition and at a relatively low speed and energy state. The 

D max 

approach is represented as: 

%C Lmax = — * 100. Equation 1 

max C L ma x 


2.2.7 Maximum Load Factor ( gumit ) 

Range for Climb: 1 .25 to 3 (g) 

Range for Descent: 0 to 0.75 (g) 

Range for Level Turn: 1.5 (g) 

These climb and descent ranges represent current normal category airplanes [7] for the 
climb maneuvers, but go to zero for descent maneuvers. Zero g is not sustainable by most 
airplanes primarily because of fuel and oil systems. For the vertical maneuvers, the 0.25 
increment on either side of 1-g straight and level flight is what TCAS expects a pilot to use for 
transition from one 1-g condition to another. Thus, it was selected as a minimum. 

The gumit effectively limits rate of change of the flight-path angle vertically (climb and 
descent) and track horizontally (level turn). This can be represented mathematically as: 

for the vertical maneuver: 


Ymax — y * id limit 1 )> 


Equation 2 
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for the horizontal maneuver: 


4>r, 


= 9 -* 


9 limit 


where g is the constant acceleration due to gravity. 


Equation 3 


Maximum load factor available may be aerodynamically limited to less than the command 
limits at low speeds due to the amount of aerodynamic lift available at the current flight 
condition, i.e., if the airplane is below its “corner speed” or “maneuver speed” it cannot generate 
enough aerodynamic lift to reach the structural g Urn i t (see [8,9]). In that case, the simulation 
models the aerodynamic gu m i t . A single level turn gu mit was used in this study. 

The load factor is expressed as: 

T —XpxV^xCj 

n z = w = - — E ■ Equation 4 

s 


The aerodynamic limitation results from a limited amount of lift an airplane can produce. 

n z max = - yiimit ■ Equation 5 


For this study, all airplanes are initially trimmed to steady level flight. At this flight condition, 
the lift the airplane generates is equal to the weight. Thus, Equation 5 becomes: 


n 7 


_ u max _ max _ 


C Ln 


Li 


CLi 


%c Lr 


— 9 limit- 


Equation 6 


Using Equation 6, in order for an airplane to generate 3 g’s, the maximum value within the 
range of interest, it must be flying at an initial C L less than or equal to 33.3% of C Lmax . 


2.3 Simplified Airplane Performance Model 

A simulation tool was developed to evaluate SAA maneuver time requirements for a 
spanning set of airplane performance models over a broad range of encounters. The tool is 
referred to as the 2-degree-of-freedom Prototyping Airplane Interaction Research Simulation 
(2PAIRS). While the 2PAIRS tool technically models more than two degrees-of-freedom of the 
airplane, the tool is limited to two-axis motion for any given maneuver, i.e., vertical climb, 
vertical descent, or level turn individually. The 2PAIRS tool was developed in the MATLAB® 
programming environment. 

The 2PAIRS simulation tool was developed to analyze the interaction between two airplanes 
involved in a single encounter. The airplane performance model is based upon simplified two- 
dimensional airfoil theory with finite-wing approximations. These assumptions allow for a simple 
model that includes performance limitations while still running fast enough for a comprehensive 
study of the overall design space. 

The simulation uses the basic design parameters from Table 1 : — , — , AR, C, . - , and 

a K W' S’ Lmax’ Dmax ’ 

Oswald’s wing efficiency factor (e). Basic aerodynamic properties are estimated from these 
parameters using two-dimensional airfoil theory [8]. This simple approach was used to provide 
generally applicable results, versus understanding the performance capabilities of a specific 
airplane configuration. Moreover, this allows evaluation of a spanning set of performance 
capabilities by varying the design parameters across a range of values of interest in a matrix 
form. To estimate airplane performance, basic airplane performance principles were used, 
starting with the assumption that straight and level flight requires thrust to equal drag (T = D) 
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and lift to equal weight (L = W). The following discussion provides some background on the 
approach. 

Aerodynamic performance metrics, i.e., lift and drag coefficients, can be estimated using 2-D 
airfoil theory applied to 3-D elliptical wings resulting in: 



r 2*n*AR*a 

L, ] — 

L AR+2 

Equation 7 

and 

P 

ii 

O 

+ 

S3 / 

Equation 8 

where a is the angle of attack in radians. The parasite drag coefficient (Q 0 ) may be 
estimated using: 


O 

II 

515 £ 
t M & 

§ S3 

Equation 9 

Then, for level flight 

//UiA 



L = qSC L = W, 

Equation 10 

and consequently 

r -- 

L ~ qs 

Equation 1 1 

where 

q = - 2 pv 2 . 

Equation 12 


The 2PAIRS simulation tool was developed using the aerodynamic performance 
approximations above and simplified equations of motion for the vertical and horizontal planes 
separately. This resulted in kinematic equations that drove the simulation as follows: 

Forces acting parallel to the flight path for a vertical maneuver: 

T cos(a) — D — W sin(y ) = m^. Equation 13 

Forces acting normal to the flight path for a vertical manevuer: 

L — Wcos(y) + Tsin(a ) = mV ^ Equation 14 

where y is the flight-path angle in degrees. Similarly for the horizontal maneuver: 

T cos(a) — D = Equation 15 

L sin(cp) = mV ^ = m^- Equation 16 

Lcos(0) = W Equation 17 


where 0 is the bank angle in degrees and 0 is the track angle in degrees. 

The high-level descriptions above are here to provide the reader insight to the modeling 
approach used and therefore fidelity of the simulation. For more information please refer to [8]. 
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2.4 Simulation Description 


The simulation initializes each run with two airplanes on a collision course, the UAS and a 
non-maneuvering intruder. The intruder may be flying level or given a constant rate of climb or 
descent. It may approach the UAS’s initial trajectory through 360° of azimuth, although a more 
limited azimuth range that accounts for right-of-way rules was used in the current study. At time 
zero, the UAS initiates the commanded maneuver and attempts the maneuver, but may be 
limited by its maneuver performance envelope. The UAS will transition to the commanded climb, 
descent, or turn rate based on gu mit or roll-rate limit. It will reach and sustain the rate as long as 
performance allows. If it has insufficient performance for the commanded maneuver, it will 
transition to its best sustainable performance over time based on its unique capabilities. For 
example, an airplane that cannot sustain a commanded climb rate will attempt to reach the rate, 
and if it reaches the commanded climb rate, it will sustain the rate as long as possible while 
losing speed. When the selected minimum speed, best-angle-of-climb speed ( V BA ) herein, is 
reached, the airplane will reduce climb rate to sustain that speed. The simulation measures the 
CPA between the UAS and intruder (Figure 2) and stops the run just after CPA is passed. The 
simulator’s graphical user interface (GUI) version provides time histories of the UAS and 
intruder states, while a batch interface compiles initial states, states at CPA, and selected other 
mid-point states of interest. The simulation uses discrete time interval integration of the related 
equations of motion at 20 Hz. 

Maneuver Initiated Performance Limited Actual CPA 

t ♦ t > 

t = 0 t = t cpa 



The encounters are set up based on initial r to the predicted CPA. Thus, the initial range is 
greater for faster-moving UAS and intruders, but the time available to avoid a collision is 
consistent between runs with the same t. The approach of having the UAS initiate its maneuver 
at time = 0 seconds makes the results independent of sensors and trackers. With this approach 
the predicted time-to-go at which a maneuver must be initiated to result in the desired CPA can 
be discerned. Regulators can use the data to help determine when a SAA system needs to act, 
define an acceptable design space, and develop required capabilities for SAA systems taking 
UAS performance into account. UAS designers can back out times for Detect, Track, Evaluate, 
Prioritize, Declare, Determine, Command and Execute (Figure 3) based on their SAA system 
design to determine when initial detection must occur. 
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Three separate 2PAIRS implementations were developed to assess basic maneuver 
elements: vertical climb, vertical descent, and level turn. The simulations account for airplane 
performance limitations associated with the particular maneuver, e.g., climb performance 
limitations in the vertical climb version or V NE in the vertical descent version. The UAS 
performance in 2PAIRS is limited by either aerodynamic/propulsive limitations, maneuver 
command limits, or a combination of both. Note that current UAS inner-loop controllers generally 
limit maneuvers far below airframe capability. The variation in commanded rates models those 
limitations in 2PAIRS as well as cases with commands greater than the UAS can sustain. 

The 2PAIRS batch capability was set up to allow for a multitude of parameter combinations 
to be run and analyzed in a time-efficient manner. This requires marching through each 
parameter variation while holding the others constant with no consideration of whether each 
particular parameter combination results in an airplane that can fly at a given initial condition, or 
even at all. Additionally, it does not make any assessment of whether an airplane would be a 
‘good’ design or not. This was deemed to be the most efficient way to set up the test matrix and 
assure that the set of results would be a spanning set of likely UAS performance within the 
region selected. This approach required some pre-run constraints to be implemented in the 
batch simulation to limit analyses to the airplanes that could fly (meeting physics-based 
constraints) and that were at an initial condition suitable for the maneuver to be executed. If 
these checks were met, the simulation would proceed with the particular parameter 
combination; otherwise, the combination would be skipped. 
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Figure 3: Diagram from 2 nd FAA SAA workshop report showing SAA self-separation and 

collision avoidance functions on a timeline. The work herein starts at “Execute.” [2] 
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The first constraint requires that the thrust available to the airplane is greater than or equal 
to that required to maintain steady, level flight at the initial knots calibrated airspeed (KCAS) 
(i.e., the initial condition of each airplane within the simulation). This constraint is represented 
mathematically by: 


T T _ T 

w max ~ w i WsiUF ' 


Equation 18 


If the thrust available is only equal to the thrust required to maintain steady, level flight, the 
airplane is able to maintain the initial flight condition, but will not have excess power to enable a 
climb without sacrificing speed. More thrust available can provide at least some climb rate 
without sacrificing speed. Conversely, if the thrust available is not at least equal to the thrust 
required, the airplane cannot maintain the straight, level, unaccelerated initial trim condition and 
will not be analyzed. 

The next physical constraint of the airplane is associated with the simulation minimum 
allowable velocity of an airplane for the desired maneuver. The constraint requires the vehicle’s 
1-g stall speed to be below V BA - Since speed is dictated by C L , this constraint may be 
represented by: 


C L ^C Ln .. 

^ max L BA 


Equation 19 


If this constraint is violated, the airplane will attempt to fly below its stall speed, which would 
result in a loss of control for the airplane. 


2.5 CPA Ranges of Interest 

Some consideration of CPA ranges of interest is required both for analyzing results and for 
scoping the test matrix. Given the test matrix approach used, many post-maneuver CPAs in the 
results would seem to far exceed the critical CPAs of interest. However, which CPA ranges are 
of interest has not yet been determined by the community. This section discusses an attempt to 
determine a reasonable range for CPAs of interest, noting that results need to include any 
possible CPA range that may be determined to be of interest in the future. The FAA 2 nd SAA 
Workshop Report [2] describes notional boundaries (Figure 4) but only associates numbers with 
the NMAC volume. The report provides descriptions but no values for the Collision Avoidance 
Threshold (CAT), Well Clear Volume (WCV), or Self-Separation Threshold (SST). Note that the 
definitions of CAT and SST include “last chance to maneuver to avoid” the next boundary. Thus, 
these two bounds are encounter and airplane-performance dependent. 

In order to draw conclusions from the results, there is a need to determine what CPA values 
would mean success. The measures of success must include both successfully avoiding a 
collision and maintaining well clear or self-separation. Thus, both small and large CPAs are of 
interest. The NMAC volume has been accepted by the aerospace community to be a cylinder 
defined as ±100 ft of the intruder vertically with a horizontal radius of 500 ft. 
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Figure 4: Figure from FAA 2 nd SAA Workshop showing collision avoidance and self-separation 

boundaries. [2] 
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The WCV is much more difficult to define. Results from this study may help define it and are 
expected to contribute directly to the definition of the CAT and SST. In order to put an upper 
bound on the goal of obtaining a spanning set of CPA results, some assumptions are required. 
The intent of these assumptions is to assure that results from this study are a superset of 
ranges ultimately determined for the WCV. That means our encounter set must include CPAs 
that exceed the currently undefined WCV. The current regulations of CFR 14 were used to gain 
some insight on what an upper limit is likely to be. CFR 14 Part 91.159 and 91.179 [1] state that 
aircraft operating under Instrument Flight Rules (IFR) on easterly headings are to fly at odd 
thousands of ft MSL, while ones on westerly headings fly at even thousands of ft MSL. Visual 
Flight Rules (VFR) aircraft are required to use the same even/odd thousands of ft altitudes, but 
add 500 ft (ee Figure 5). The current airspace allows these aircraft to pass with these vertical 
separation distances and no horizontal separation. Thus, it seems reasonable to assume that 
either 500 or 1000 ft vertical separation is the current definition of vertical well clear. As a result, 
this study presents results in relation to three values of vertical separation: 1 ) < 1 00 ft for NMAC; 
2) < 500 ft referred to as IFR-VFR; and 3) < 1000 ft referred to as IFR-IFR or VFR-VFR. CPAs 
greater than 1000 ft are considered to have exceeded the WCV requirements likely definition. 
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Figure 5: Information used 

A maximum horizontal CPA is more difficult to discern. Controllers are required to separate 
IFR traffic by 3 nautical miles in terminal areas. However, VFR traffic self separates visually and 
even two aircraft operating under IFR may be allowed to self separate after the pilots 
acknowledge to controllers that they have visually acquired the other aircraft. The result is that 
self separation is maintaining well clear as defined by the onboard pilot’s discretion. Thus, well 
clear could mean anything inside 3 nm. As a result, this study will assess horizontal separations 
encompassing 500 ft (NMAC) out to 3 nm (terminal area IFR separation). 


Thus proposed spanning set of separation distances 
covering NMAC through Self Separation Threshold is 
expected to be no greater than: 

1 00 feet < CPA < 1 000 feet vertical separation 
500 feet < CPA< 3 nm horizontal separation 

to Select Maximum CPAs of Interest. 
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2.6 Avoidance Maneuvers and Encounter Parameters 


This section describes the avoidance maneuvers and the intruder behavior modeled. 
Parameter ranges are given for the commanded vertical climb and descent and horizontal turn 
avoidance maneuvers, and for the encounter scenario and intruder behavior. Descriptions of the 
GUI simulator tools are also given. 

2.6.1 Climb Maneuver 


When a simulation run starts, the UAS is trimmed in straight and level flight and immediately 
initiates a maneuver to the commanded climb rate. The UAS transitions from level flight to the 
commanded climb rate at the maximum load factor (referred to as max nz in Figure 6) set for 
the maneuver or the aerodynamic limit, whichever is lower. If the UAS has sufficient 
performance to sustain the climb without decelerating, it will do so. If it does not, it will lose 
speed while maintaining the climb rate (if possible) until the minimum allowable speed, or V BA , is 
reached. At that point it will push over to the corresponding sustainable climb angle using the 
same g-increment used in the pull up. It will then continue the climb at the sustainable climb 
speed and angle. 

The GUI for the vertical climb maneuvers is shown in Figure 6. The user can enter numbers 
in each box to visually see the results of changes. By clicking the button next to each parameter 
name, the slider can be used to see the resulting trajectory and CPA change as the value of the 
selected parameter changes incrementally. Results can be sent to the workspace to allow 
analysis of time histories of a wide number of states for the run. Additionally, the colored 
squares allow the plot to be exported. Selecting another color after changing parameters will 
cause the result to be plotted with the result from the last time a colored button was selected. 

Up to six results can be added to a plot in this manner. The batch tool allows sweeps of all 
parameters shown on the GUI and outputs initial, final, and some interim states of interest for 
the UAS and intruder. 
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Figure 6: 2PAIRS GUI tool - Vertical Climb. 
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2.6.2 Descent Maneuver 


Similar to the climb maneuver, the UAS is initially trimmed in straight and level flight. When a 
simulation run starts, the airplane immediately initiates the commanded descent rate. It rotates 
at the commanded or minimum aerodynamic load factor for the configuration to reach the 
commanded descent rate. The airplane attempts to maintain the initial speed by reducing the 

throttle position until minimum power is reached. In 2PAIRS, that means ^ = throttle = 0. 

Depending on the magnitude of the commanded descent rate, the airplane may gain speed 
while continually adjusting flight-path angle to maintain the steady rate of descent until Vne is 
reached. If Vne is reached, the airplane will pull up at the commanded or available gu m i t to the 
glide path that will maintain Vne- The airplane will maintain this speed and descent rate until 
CPA is reached and the run ends. The GUI for the vertical descent maneuver is shown in Figure 
7. A batch tool front end was also developed. 
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Figure 7: 2PAIRS GUI tool - Vertical Descent. 
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2.6.3 Vertical Maneuver Parameters 

This section gives ranges and descriptions of parameters associated with the vertical climb 
and descent maneuvers. 

2.6.3. 1 Commanded Vertical Rate (h command ) 

Range: 500 to 5000 (ft/min) 

This is the commanded rate of climb or descent, for the respective maneuvers, that the 
airplane receives at time = 0 seconds in the simulation. The airplane will attempt the 
commanded rate while trying to maintain forward speed. If forward speed cannot be maintained, 
the airplane will adjust its flight path to maintain the limiting speed, as described earlier. 
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2.6.3. 2 Best-Angle-of-Climb Speed (V BA ) 


Range: Varies based on other parameter values 


Unique to the vertical climb maneuver, this parameter is used to restrict the minimum speed 
for the maneuvering airplane to the speed associated with - . This restriction keeps the 

D max 

airplane at a best sustained climb condition and prevents stall. If the airplane must sacrifice 
forward speed to achieve or maintain vertical speed, it will do so until V BA is reached. The 
airplane will then adjust the climb rate to the sustainable best-angle-of-climb rate associated 
with the airplane configuration. 

Vba is calculated from airplane design parameters at the start of each run: 


Vba = 


_ W L 
2 * — * — 

S Dmax 

p*n*e*AR 


* sin 


-il L 

w 



Equation 20 


In real world practice, pilots rarely fly below V BA , often referred to as loiter speed, unless 
practicing slow flight procedures and stall recovery. 


2. 6. 3. 3 Velocity Not to Exceed (V NE ) 


Range: 1.15 * V cruise (KCAS) 


Unique to the vertical descent maneuver, this parameter is used to define a maximum 
velocity the airplane can reach. This can be a structural limit, a control law limit, or a velocity 
restriction imposed by rules of the NAS (i.e., airplane speed limit). This value is not directly 
related to aerodynamic and propulsive performance, so some gross assumptions were made 
about margins above cruise speed. The results presented herein are from one increment above 
cruise. For this study, cruise speed is defined as the 1-g speed associated with 8% of the C Lmax , 
expressed as: 


V„ 


2 * 


w 


p*.08*C Lr 


Equation 21 


2.6.4 Level Turn Maneuver 

Similar to the vertical maneuvers, the airplane is initially trimmed to straight and level flight. 
At the start of the simulation run, the airplane begins to roll at the specified rate to the bank 
angle that will yield the commanded turn rate. If the commanded rate is greater than the rate of 
turn achievable within the structural limit or bank-angle limit, the airplane will be limited to the 
maximum achievable turn rate. The UAS will continue the turn until CPA is reached or the point 
at which it rolls to wings-level at the specified roll rate to achieve the commanded track change 
of 90°. It will then continue to fly level at a steady speed until CPA is reached and the simulation 
ends. The GUI for the level turn maneuver is shown in Figure 8. A batch front end was also 
developed for this maneuver. 
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Figure 8: 2PAIRS GUI tool - Level Turn. 
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2.6.5 Level Turn Parameters 

This section gives ranges and descriptions of parameters associated with the level turn 
maneuver. 

2. 6. 5. 1 Bank Angle Limit ((t>i imit ) 

Range: 10° to 60° 

This parameter limits the bank angle an airplane can achieve. This parameter has a direct 
effect on an airplane’s achievable turn rate. There is a direct relationship between the bank limit 
and imposed gu m i t that is not related to airplane speed. However, the bank angle or gii mit 
determine the turn rate and radius for a given airplane speed. For this study, gu m i t is used to set 
the bank angle limit. The equation relating these two parameters can be expressed as: 

Sumit = E « ua,ion 22 

2.6. 5.2 Heading Change Command (Atp) 

Range: 90° 

This parameter is the commanded heading change the airplane will make before returning 
to straight and level flight. The implementation makes this heading capture within ±2° for all 
runs. 
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2.6. 5.3 Roll Rate Limit ((/) limit ) 

Range: 30 (deg/sec) 

This parameter sets the roll rate for transition between straight and level flight and the level 
turn and back to straight and level flight. This value is not expected to be a major player given 
the required time associated with horizontal turns, i.e., the time required to roll in and roll out of 
the turns is small compared to the time to turn to the specified heading. Thus, only one value 
was run in the dataset. Additionally, the airplane reaches and stops the roll rate in one time 
step, i.e., roll acceleration is not modeled. 

2.6. 5.4 Turn Rate Command (iff command) 

Range: 1 to 10 (deg/sec) 

This parameter is the commanded turn rate the airplane will attempt to maintain. The 
airplane may not reach the commanded rate if there are turn rate limits imposed by a bank limit 
or gu mit , as outlined in Equation 3. 

2.6.6 Encounter Parameters 

This section gives ranges and descriptions of the encounter parameters used with all three 
maneuvers. The encounter geometry definitions, azimuth in particular, align with a white paper 
presented and discussed in RTCA SC-203 and that will be discussed in SC-228. It has not yet 
been formally published. However, the definitions are described in [10]. 

2. 6. 6 . 1 Time-to-go (x) 

Range: 10 to 120 (seconds) 

Time-to-go is the predicted time until CPA between two airplanes involved in an encounter. 
For the encounters simulated, the predicted CPA was always 0 ft, i.e., a collision. If no 
maneuver is performed, the two airplanes will collide or reach CPA x seconds after the run 
starts. This parameter is also used to define the initial range between the two airplanes given 
the initial velocity of both airplanes and the encounter geometry. It should be noted that the 
initial time-to-go (predicted x) is not necessarily the time at which the actual CPA occurs due to 
changes in the closure rate and encounter geometry as a result of the UAS’s maneuver. Since 
the initial x range must encompass Collision Avoidance and Self-Separation thresholds, x was 
varied from 10-120 seconds in the simulation. 

2. 6. 6. 2 Intruder Velocity (V intr ) 

Range: 50 to 250 (KCAS) 

The intruder flies at a constant velocity throughout the simulation, input as KCAS. The 
speed of intruders ranges from 50 to 250 KCAS for the full range of altitudes tested. This is the 
speed limit in the NAS below 10,000 ft MSL. Thus, faster intruders could be encountered above 
that altitude. Aircraft above 10,000 ft MSL are required to have transponders so they are very 
likely to be dealt with by ATC. Still, faster intruders should be assessed at a later date. 
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2.6.6.3 Intruder Azimuth (fi) 

Range for Vertical Maneuvers: 0° to 90° 

Range for Horizontal Maneuvers: -90° to 90° with UAS always turning right. 

The angle of the encounter was limited to the forward quadrants of the UAS in order to avoid 
scenarios with the UAS being overtaken. These scenarios were omitted because, in the current 
airspace, it is always the duty of the overtaking aircraft to avoid the slower-moving aircraft, and 
this rule is currently expected to carry over to UAS. The intruder was always heading toward the 
UAS as well, thus overtaking scenarios by either the intruder or the UAS were not included in 
the results to date. They would require lower maneuver performance than the scenarios 
evaluated. 

From Figure 9 and reference [10], /? is the angle measured between the ownship’s and the 
intruder’s initial horizontal flight paths. For vertical maneuvers (i.e., vertical climb or descent), 
this value will range from 0° (head-on) to 90° (crossing) because results will be symmetrical for 
approaches from the other side. The horizontal maneuver ranges from -90° (crossing from left) 
through head-on to 90° with the UAS always turning right. Symmetry can be used for assessing 
left turns. 

For this study, the range at CPA ( R CPA in Figure 9) is always 0 ft, meaning the two airplanes 
always begin on a direct collision course. Thus, there are no asymmetries associated with 
encounter offsets. 


2.6. 6. 4 Intruder Vertical Rate (h intr ) 




r c 




R , = V, * r„ 



Range: -5000 to 5000 (ft/min) 

Intruder vertical rate is input as a positive or 
negative rate in ft/min. The value is held constant 
throughout a simulation run. It is assumed that the 
intruder can maintain the specified vertical rate, 
and therefore, no performance limitations are 
considered for the intruder. The vertical rates 
simulated for the intruder range from a constant 
climb at 5000 ft/min to a constant descent at 5000 
ft/min. While the full range of intruder climb and 
descent rates were run, a post-analysis exclusion 
was placed on the maximum flight-path angle from 
which the intruder could approach. For the data 
presented herein, the intruder did not approach at 
a vertical flight-path angle with magnitude greater 
than ±12°. For instance, this limitation excludes 50 
KCAS intruders climbing or descending at 5000 
ft/min, which results in an extreme and 
unreasonable flight-path angle of 80°. 


Figure 9: Encounter Geometry and 
encounter parameter 
definitions. Diagram from 
SC-203 document. [9] 
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2.6.6. 5 Altitude (h) 


Range: 1000 to 18,000 (ft MSL) 

Measured in ft MSL, altitude is used primarily to specify the air density for the encounter. 
During the run, avoidance maneuver vertical deviations of interest are small. Thus, atmospheric 
density is held constant through each simulation run to save computational complexity. Test 
runs were conducted to show this assumption made no significant difference in the resulting 
trajectory or CPA. Altitude has a direct effect on the difference between the calibrated airspeed, 
which drives aero-propulsive performance, and true airspeed, which dictates the inertial 
reference frame. The test range chosen for the encounters ranges from 1 ,000 to 18,000 ft MSL, 
which is generally landing pattern altitude to the bottom of Class A airspace. This range follows 
the RTCA SC-228 Terms of Reference which outlines that the focus of initial work will be on 
UAS transitioning to and from Class A airspace [4] or special purpose areas. 

2.7 Verification and Validation of Simulation 

Before moving into the large volume of data available from the 2PAIRS batch tool, data from 
individual airplanes found in Table 1 were used for verification and validation of results from the 
simulator and to illustrate trends. Figure 10 shows climb performance results for a 
representation of a Cessna 172, which is a single-engine general-aviation airplane that has four 
seats and is capable of cruising at approximately 110 KTAS and achieving a best rate of climb 
of approximately 750 ft/min at maximum gross weight and MSL. Results shown are from runs 
conducted at 1000 ft MSL with x = 30 seconds using the GUI tool which provides time history 
data. The simulation data for the Cessna 172 matched the behaviors that were anticipated, 
showing that the simulation accurately modeled how the Cessna airplane could react in the 
scenarios. At time zero, the airplane was commanded to climb at 5000 ft/min, which is well 
beyond its capability. The four runs represent the airplane at near maximum gross weight (high 

y and low and near a minimum likely flying weight (low j and high ^). The low-speed initial 

conditions are at Vba speed and the higher-speed conditions represent cruise speed. The 
results show that the airplane immediately went to Vba climb condition from the two slower initial 
conditions. For the faster initial conditions, the airplane tried to achieve 5000 ft/min but was 
performance limited as expected. As it lost speed, the climb became steeper until Vba was 
reached. Then climb rate was reduced to the rate of the sustainable condition. Note the CPA 
differences. The worst-case, heavy and slow configuration, took approximately 9 seconds to 
clear an NMAC and did not achieve the IFR-VFR separation of 500 ft in 30 seconds. The other 
three configurations resulted in clearing the 500-ft separation threshold. Only the faster and 
lighter initial condition exceeded the 1000-ft CPA threshold. 
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1000 
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Figure 10: Data from Cessna 172 using 2PAIRS climb tool. 


Figure 11 shows CPA results from all of the configurations in Table 1 as a function of y, 

Vi, and specific excess power ( P s ). Results from each configuration with both a heavier and 
lighter weight condition are shown. Time histories of these scenarios were inspected as part of 
verification and validation but are not presented in this paper. Note that no attempt was made to 
model the flight controls of any of the aircraft in this work, only the aerodynamic/propulsive 
potential to maneuver was modeled. Additionally, neither sensor noise nor navigation error were 
modeled; thus, the results are for idealized navigation. 

The largest CPA values in Figure 1 1 create a fairly flat line across all subplots. This is a 
result of the same commanded climb rate being used for all runs and is representative of the 
configurations that could achieve and nearly sustain the rate. With insufficient performance to 
achieve or sustain the commanded climb rate, the results varied. The resulting CPAs for this 
insufficient performance group show no definite trend against parameters other than P s , which is 
a key performance parameter related to vertical maneuvers. The scales are not shown on the 
plots because the reader should not assume the results represent the actual aircraft, especially 
the UAS. In the simulations presented, a high commanded climb rate was input and the 
modeled aircraft tried to follow it. Most current UAS control laws would not allow such a high 
rate to be attempted. Also, while the trends look very reasonable, the simulation is not intended 
to match a specific airplane but rather provide a broad span of performance results. 
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Figure 1 1 : CPA Results For Heavy and Lightweight Conditions for All Configurations in 
Table 1 as a Function of y, V h and P Sr Scales are not shown on the plots 

because the results should not be assumed to represent the actual aircraft, 
especially the UAS. 


3 Results and Analysis of Climb Maneuver 

This section presents the data from the vertical climb maneuvers along with analysis of the 
results. Descriptions of key influential performance parameters are included as they are 
presented. Data from the vertical descent maneuvers is presented in Section 4, and level turn 
maneuver results are shown in Section 0. 

3.1 Climb Maneuver Performance Parameters 

Analysis of this multi-dimensional design space is an iterative process with a goal of 
characterizing the relationship between design and performance parameters to CPA values 
across the range of airplanes and encounters. In order to deliver the most usable information for 
both UAS designers and the regulators developing certification requirements to allow UAS 
regular access to the NAS, performance parameters that are readily available, practical, and 
meaningful were sought. Maneuver performance for a given encounter is a function of the 
airplane’s design parameters and initial condition. Thus, performance parameters which are 
derived from the design parameters and initial flight condition were anticipated to be more useful 
than design parameters alone because they relate more directly to the airplane’s 
maneuverability and, therefore, its ability to increase CPA for a given encounter. 

For reference during the following discussions, the design and encounter parameters 
relevant to the climb maneuver simulations are listed in Table 2 showing the discrete values 
simulated. These parameters were described above in Sections 0 and 0. 
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Table 2: Simulation Parameters Relevant to Climb Maneuver 


Parameter Name 

Parameter Symbol 

Discrete Values Simulated 

Wing loading (lbs/ft 2 ) 

w 

~s 

1,9, 17, 25 

Thrust-to-weight ratio 

T 

w 

0.1, 0.2, 0.3, 0.4 

Aspect ratio 

AR 

5, 16.67, 28.3, 40 

Maximum coefficient of lift 

^ L max 

0.6, 1.3, 2 

Initial percent of maximum 
coefficient of lift 

%C L maX 

8, 44, 80 

Maximum lift-to-drag ratio 

L 

D max 

5, 20, 35, 50 

Load factor Limit 

9 limit 

1.25, 3 

Altitude (ft MSL) 

h 

1000, 6666.7, 12333, 18000 

Commanded climb rate (ft/min) 

h command 

500, 2000, 3500, 5000 

Intruder velocity (KCAS) 

Vintr 

50, 150, 250 

Intruder vertical rate 

hintr 

-5000, -2750, -500, 0, 500 
5000 

Intruder azimuth (deg) 

P 

0, 45, 90 

Time-to-go (seconds) 

T 

10, 20, 30, ... 120 


3.1 .1 Specific Excess Power ( P s ) 

The first step was to consider the total set of climb maneuver data to look for overall trends 
and start to drill down to useful results. Based on previous experience and the verification and 
validation study data shown in Figure 1 1 it was anticipated that P s , or some related parameter, 
would be a good predictor of climb maneuver performance. The energy state of an airplane is 
described by the combination of speed and altitude: kinetic and potential energy, respectively. 
Specific excess power is a measure of the ability to change energy state [9], i.e., increase or 
decrease velocity, altitude, or a combination of both. Thus, the performance parameter P s 
reflects the maneuverability margin of a given airplane from a given flight condition. Specific 
excess power can be represented mathematically in many forms. A relatively simple one was 
suitable for this work: 


P s = *V = V sin(y) = h. Equation 23 

This equation can be expanded to directly use the design parameters outlined in Section 0: 


p _ I T ireAR Cl 

* — \W .I 2 r neAR 

4 d Cl 


il 

I P C L 


Equation 24 


From Equation 23, it is clear that if P s is a positive quantity, that is to say thrust is greater 
than drag, the airplane is able to increase energy, i.e., climb at some rate without losing speed 
(change in potential energy) and/or accelerate in level flight (change in kinetic energy). Note that 
an airplane that is at or near maximum throttle and speed can only climb by sacrificing kinetic 
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Figure 12: Design Space Overview; Commanded Climb Rate Shown by Color. 


energy for potential energy and would have a P s value near zero. Also, for an airplane trimmed 
to straight, level, unaccelerated flight, the excess power available to change from the initial 
energy state is zero if the throttle is not changed. The airplane must adjust throttle position to 
either accelerate or climb without losing speed, and thus, the initial P s is not as useful as the 
potential of the airplane to change energy state. For this reason, P s calculations are typically 

made using maximum ^ available. Under this assumption, the equation for the initial P s for the 
vertical climb becomes: 


Pc 


( L 

\Wmax 


neAR 


Cl 

ne AR 


2 I 

P Cl 


Equation 25 


For the discussion of vertical climb performance, all references to Pc. refer to the Pc. , as 
calculated in Equation 25. 

Figure 12 shows the relationship between P s . and CPA for all 23.27 million climb maneuver 
runs resulting from the matrix of data in Table 2. A series of horizontal lines that represent a 
common CPA response regardless of the P s . value can be seen in the figure. These horizontal 
lines are from multiple UAS configurations that have sufficient P s . to maintain the commanded 
climb rate without sacrificing forward speed. As indicated by the color legend, lines correspond 
to a constant h command and x.This is discussed in more detail in the performance group 
discussion that follows. 

Clear vertical lines are also apparent and become further separated as P s increases. These 
vertical lines are grouped by close P s . values. For instance, there is a pair of vertical lines at 
constant P s . values of -138 ft/s and -140 ft/s. The UAS represented by these P s . values have 
common design parameter values except for - . The higher P s . value relates to a higher - 

D max 1 D max 

value because as - increases, the UAS requires less of the available power to maintain 

Dmax 

straight, level, unaccelerated flight. It therefore has more excess power at the initial trimmed 
condition. 
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Figure 13 shows the experiment design space separated by r. It is worth noting that there 
are bands of horizontal lines that share an h command (i.e., the horizontal lines are the same 
color). The CPA difference between the top and bottom of the bands are due to variations in 
gumit > ar| d are more apparent than in the previous figure. 


CPA (ft) vs. Specific Excess Power (ft/s) 
Tau(sec) 

40 50 60 70 



'88 ° 8 § ° 8 § ° 8 § 
Specific Excess Power (ft/s) 


8 ° 8 8 ° 8 8 ° 8 S 


Commanded Climb Rate (fpm) 

500.0 

1000.0 

1500.0 

2000.0 

2500.0 

3000.0 

3500.0 

4000.0 

4500.0 

5000.0 


Figure 13: Design Space Overview; Segregated by Time-to-go; Commanded Climb Rate 
Shown by Color. 


3.1 .2 Analysis of Specific Excess Power’s Effect at Various Flight Conditions 


One relationship of interest is the difference between the airplane configuration’s P s . and 
specific excess power at the Best Angle of Climb, ^ 'sba ■ This relationship can be expressed 
mathematically as: 


^ ' S BA — ^ s i BA — h Sustainable h BA ■ 

Equation 26 

where h Sustainame is the maximum sustainable climb rate at the airplane’s initial speed. The 
relationship can be conditionally simplified to: 

If ba ^ then hs us tainable ^ h BA 

Equation 27 

If AF 'sba ~ then fysustainafr/e = h ba 

Equation 28 

If AF *Sba ^ then hs us t a inable ^ h*BA ■ 

Equation 29 


After studying these and the underlying relationships, it can be shown that it is not valid to 
exclude any range of A P Sba data from the study. Since the Best Angle of Climb is an 
optimization of the flight-path angle rather than the rate of climb, there is no constraint on 
h sustainable with respect to h BA . Therefore, A P Sba effectively represents the relationship between 
h. sustainable an d h BA . As a result, there is no constraint on A P Sba . 
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3.1.3 Potential Change in Velocity (AV) 

For the climb maneuver, this parameter is a direct measure of the speed (kinetic energy) a 
given airplane can lose from the initial condition before reaching the limiting velocity associated 
with the selected maneuver, i.e. V BA . It is essentially a measure of the kinetic energy that can be 
traded for potential energy (altitude) before reaching the performance-limited state. It is 
calculated as: 


AV = V t - V min . Equation 30 

As the initial velocity of the airplane approaches the minimum velocity, the available speed 
margin decreases to zero. This affects the amount of time an airplane has to achieve, or that it 
can maintain, the commanded climb rate before adjusting to a sustainable climb condition. As 
the value of AV increases, the amount of time the aircraft can maintain the commanded rate 
before adjusting to a lesser rate increases for common deceleration rates. More time at the 
higher climb rate increases the achievable separation during an encounter. 

3.1.4 Difference between Commanded Climb Rate and Sustainable Climb Rate 

(Ahsustainable) 


This parameter reflects the rate of speed lost by the airplane for the commanded climb rate. 
The sustainable climb rate at constant speed of a given airplane is directly measured by P s .. The 
equation that defines this parameter is: 

Sustainable = h command ~ ^S, ■ Equation 31 

If the commanded climb rate is less than the sustainable rate, this value is negative, and the 
airplane is able to achieve and maintain the commanded rate without sacrificing forward speed, 
which results in geometrically predictable CPA values. As this value approaches zero (i.e., 
commanded climb rate approaches the sustainable climb rate), more throttle is required to 
maintain the commanded rate at the initial speed. If this value is positive, the commanded rate 
is greater than the airplane’s sustainable climb rate. This means the airplane will be forced to 
sacrifice speed to maintain the commanded rate. As the simulation continues, the airplane’s 
speed will approach V BA and it may have to adjust to the associated sustainable climb rate. The 
higher the value of Ah Sustainablet the faster the aircraft will decelerate to maintain the 
commanded rate. Combined with AV outlined in Section 3.1.3, the magnitude of this parameter 
has a large effect on achievable CPA. 

3.2 Performance Grouping of Climb Maneuver Data 

Numerous attempts to develop response surface models from the entire data set did not 
result in sufficiently accurate models. Thus, the data was divided into more homogeneous 
groups in a search for more useful alignment of analysis. Several different ways of grouping the 
climb maneuver data were considered. Initial response surface models were fit to preliminary 
data and residual analysis was done. Results contained several groups of responses 
independent of changes in airplane parameters. These groupings were revealed to be 
negatively affecting the accuracy of the model fit. Using JMP® [11] data exploration tools and 
2PAIRS GUI functionality, individual airplane configurations within each grouping were 
analyzed. 

Four major groups were developed, segregating airplanes by climb maneuver performance 
capabilities and resulting in more accurate model fitting for the airplane of interest. Notional 
climb maneuver performance representative of each of the four groups is shown in Figure 14. 
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The Sufficient Power Differential (SPD) group represents cases where airplanes have sufficient 
power to maintain the commanded climb rate throughout the given encounter without sacrificing 
forward speed. The Insufficient Power Differential (IPD) group is where airplanes do not have 
sufficient power to maintain the commanded climb rate without sacrificing forward speed. In 
these cases, the airplanes initiate the commanded climb, but as their speed slows, they must 
reduce their climb rate in order to not go below V BA . In the Insufficient Initial Speed (IIS) group, 
the airplanes have insufficient initial speed to conduct the commanded climb, i.e. their initial 
forward velocity is below V B a- These airplanes hold a level flight path and increase power output 
until V BA is reached and then proceed to climb at the steady climb rate associated with V BA . The 
Insufficient Time to Pitch (ITP) group represents cases where airplanes have insufficient time to 
pitch up to the commanded climb rate before the CPA is reached, thus they are still in the 
process of pitching up to the commanded climb rate throughout the entire encounter. Figure 15 
shows CPA vs. P s . for each of the performance groups and highlights the differences between 
these groups. The climb maneuver performance for each of these groups is considered 
separately in Sections 3.2.1 through 3.2.4. 



Figure 14: Notional chart depicting four climb maneuver performance groups 

In order for the results of the study to be generally useful, it was important to select 
performance groups that support accurate prediction of the group to which a particular airplane 
configuration and flight condition belongs. 2PAIRS has the ability to predict the performance 
group of an airplane configuration before running the time history. An airplane configuration 
alone does not determine the performance group. The combination of airplane parameters, 
commanded maneuver rate, and initial flight condition are all required to determine performance 
group as it relates to the SAA system requirements. Additionally, small deviations within an 
airplane’s design parameters or the commanded maneuver may result in that airplane changing 
performance groups. For instance, for a given airplane combination belonging to the SPD group 
(Section 3.2.3), an adjustment in commanded climb rate above the configuration’s sustainable 
rate would result in that airplane shifting to the IPD group (Section 3.2.4). This can be 
addressed by studying an airplane configuration at different stages of its flight envelope 
performing a variety of maneuvers. 
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Each group is discussed separately in the following sections. Two of the groups are 
discussed in detail including analysis of results. One of these groups produces succinct CPA 
results as a function of a few parameters. The other proves much more difficult but some 
conclusions are drawn and a methodology is presented. The other two groups are explored in 
far less detail because it is unlikely either represents a real-world case. The two that are of less 
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Figure 15: Design Space Overview; Segregated by Group; Time-to-go Shown by Color. 


interest are explored briefly and first. It is worth reinforcing that the airplane design parameters 
and resulting performance parameters alone do not define the group. Thus, most design 
parameters that represent an airplane that can fly at all are represented in some part of the two 
groups that are explored in detail. 

3.2.1 Insufficient Initial Speed Group 

This airplane performance group is defined by airplanes whose initial velocity is below V BA 
as described in Section 2. 6. 3. 2. The simulation logic addresses this situation by forcing the 
airplane to hold a level flight path while accelerating to V BA . The airplane then proceeds to climb 
at the steady climb rate associated with V BA - These airplanes were not studied in detail because 
it is unlikely that an aircraft would fly in this condition under normal circumstances. Additionally, 
the same airplane configurations that occupy this performance group will also be given higher 
initial velocities as other test matrix points, which will place them into the performance 
categories of interest. 

3.2.2 Insufficient Time to Pitch Group 

This performance group is characterized by airplanes that are unable to reach the 
commanded climb rate due to very short simulation times and/or very high speeds. This group 
can be modeled using simple geometry and known performance limitations (i.e., limited pitch-up 
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rate). Note from Figure 15 that this performance group contains a very small number of 
simulation configurations. Specifically, it only includes encounters with t of 10 seconds. 

3.2.3 Sufficient Power Differential Group 

Airplanes able to maintain the commanded climb rate throughout the encounter without 
sacrificing forward speed fall within this performance group. Considering the performance 
limitations of the pitch-up maneuver, it is easy to predict the trajectory of these airplanes using 
geometry. A parametric equation has been developed using response surface methodology that 
is able to accurately predict the achievable CPA for airplanes within this group. 

The airplane’s ability to maintain the commanded climb is dictated by the airplane’s specific 
excess power. As shown in Section 3.1.1, the P s . of an airplane from a particular flight state 
indicates the airplane’s ability to change energy state. Climbing at constant speed, acceleration 
or deceleration in level flight, and descent at constant speed are examples. If an airplane’s P s . is 
greater than the commanded climb rate, the airplane will be able to sustain that rate without 
sacrificing forward velocity (i.e., kinetic energy). Therefore, this performance group is comprised 
of all airplane combinations with P s . greater than or equal to the commanded climb rate. 

Since all airplanes in this group are able to maintain the commanded maneuver, there is no 
difference in achievable separation for airplanes with common design parameters but different 
thrust capabilities. For instance, an airplane that is able to sustain the commanded maneuver 

with a ^ = 0.2 will have the same result as an airplane with the same design parameters but a 
= 0.4. The more powerful airplane will not require as much throttle to sustain the maneuver, 
however. The higher ^ airplane will be capable of flying a more extreme maneuver than the low- 

power airplane so the h command becomes the limiting factor and equalizer. As a result, 
commanded climb rate has a large impact on the distribution of airplanes within this group. More 
airplane configurations in the data set, as in the real world, are able to sustain the 500 ft/min 
climb than are able to sustain the 5000 ft/min climb. Conversely, all airplane configurations 
capable of flying at 5000 ft/min are able to sustain a 500 ft/min climb. This results in an increase 
in the number of airplane configurations in the SPD group as the commanded climb rate 
decreases. 

Figure 16 shows the CPA response versus the initial range to the intruder for airplanes 
within this performance group. The colors represent the commanded climb rate. There are 4 
commanded climb rates in the dataset as shown in Table 2 and evident in the figure. This figure 
shows that, as expected, an increase in commanded climb rate results in an increase in CPA for 
a given t. 
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Figure 16: SPD Group; Commanded Climb Rate Shown by Color. 


3. 2.3.1 Analysis of Parameter Effects on Climb Maneuver - SPD Group 

Although the performance of the UAS within this performance group is highly predictable, 
various components of the encounter have an effect on achievable separation or CPA. For 
instance, the speed and direction from which the intruder approaches have an effect as well as 
the altitude of the encounter. Analysis of a number of these parameters is presented in the 
following subsections to provide a better understanding of the effects of each parameter and 
their implications on the predictive equation. 

3. 2.3. 1.1 Effect of g limit 

To show the impact an airplane configuration’s gu m i t has on achievable separation, Figure 
17 shows airplanes within the SPD group involved in an encounter with r = 60 seconds. For 
each commanded climb rate in the figure, two near-constant response values are seen due to 
the two discrete gu m i t values in the data, as detailed in Section 2.2.7. The gu m i t effect appears 
to increase as the commanded climb rate increases. For an airplane flying at a fixed speed, a 
higher commanded climb rate results in a higher commanded flight-path angle. As the 
commanded climb rate increases, the time required to pitch at a constant rate to the 
commanded rate or angle increases. With an increase in pitching time, the effect of gu m i t 
becomes more influential on the achievable separation. Upon applying Equation 2, Table 3 
shows the time spent pitching to 500 and 
5000 ft/min at several gumit values for an 
aircraft flying at 100 knots. More than 7 
additional seconds are needed for the 
airplane with limited pitching capability 
(giimit °f 1 -25) than are required by the 
more agile airplane to reach 5000 ft/min 
climb rate. Within those 7 seconds, the 
more agile airplane can gain over 550 ft. 


Table 3: Time to pitch at various gumit values for an 
aircraft at 100 knots. 


9lim.it 

Time to 
reach 500 
fpm 

Time to 
reach 5000 
fpm 

1.25 

1.03 s 

8.5 s 

2 

0.26 s 

2.13s 

3 

0.13s 

1.06 s 
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Figure 17: SPD Group; r= 60 seconds; Commanded Climb Rate Shown by Color 


Ignoring time to transition from level flight to the steady climb, an airplane climbing at 5000 
ft/min for 60 seconds is able to gain 5000 ft of altitude. During the simulation, the UAS is initially 
trimmed to level flight and given a commanded climb rate. The faster the airplane is able to 
maneuver from level to the commanded climb rate, i.e., the higher the gu m i t , the closer to the 
ideal change in altitude it can achieve. This is evident in Figure 18 as the higher gi imit value of 
each commanded climb rate is approaching the respective expected value. 



Figure 18: SPD Group; t = 60 seconds; g Um it Shown by Color. 
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3. 2. 3. 1.2 Effect of Intruder Velocity 

For a given UAS configuration, the intruder velocity in combination with r defines the initial 
range between the intruder and the point of collision. To identify the effect of the intruder’s 
velocity, Figure 19 was developed. This figure shows airplanes in the SPD group that were 
commanded to climb at 500 ft/min. Each horizontal line represents a r value from 10 to 120 
seconds. The color in the figure represents the intruder’s constant velocity. Trivially, as the 
intruder’s velocity increases the initial range between the two airplanes increases, as shown in 
Figure 19. Conversely, the initial range has no impact on CPA. The CPA is a function of r only 
in this figure. Thus, if the UAS initiates the same commanded climb rate to avoid either a fast or 
slow intruder at the same r, the achievable separation is the same despite the larger range at 
maneuver initiation. 
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19: SPD Group; Commanded Climb Rate of 500 ft/min; Intruder Velocity Shown by 
Color. 


From Figure 19, it can be seen that for a given CPA, there exists large variability in initial 
range for intruders of common velocity. Initial range is measured as a slant range between the 
two airplanes; the variability is due to the various approach azimuths, elevation angles, and 
closure rates between the two airplanes. There is also a bit of variability in CPA in the lower 
initial range values for each intruder velocity. The next subsection identifies the causes of this 
variability, focusing on a single r value. 

3. 2. 3. 1.3 Effect of Intruder Climb Rate 

Figure 20 focuses on a single commanded climb rate of 500 ft/min and r of 60 seconds. The 
diimit is constant at 1 .25 g’s for the UAS. The color in Figure 20 shows the vertical rate of the 
intruder. The lines defining ±5000 ft/min begin at a range of about 25,000 ft MSL due to the 
limitations in intruder climb angle described in Section 2. 6. 6.4. From the figure, it appears that 
the magnitude of the intruder’s vertical rate is more significant than the sign. This is evidenced 
by the pairing of vertical rates with common magnitudes. It is important to note that the 
variability in response in Figure 20 is about 30 ft, 6% of the predicted separation. That difference 
can likely be deemed insignificant in practice. 
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Figure 20: SPD Group; r = 60 seconds; gu m i t = 1.25, Commanded Climb Rate = 500 ft/min; 
Intruder Vertical Rate Shown by Color 


To identify variance among encounters with common intruder vertical rates, Figure 21 shows 
only the ±2750 ft/min rates for the same configurations shown in Figure 20. From the figure, it 
appears that the descending intruder creates more separation than the climbing intruder. There 
are also clear groupings of performance trends for each vertical rate. Again the variation is 
small. Both of these were explored to understand fingers seen in the larger dataset. 



Figure 21: SPD Group; r = 60 seconds; gu m i t = 1.25 g’s; Commanded Climb Rate = 500 ft/min; 
Intruder Vertical Rate = ±2750. 
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3. 2. 3. 1.4 Effect of Intruder Azimuth 

Figure 22 shows the CPA response versus initial range for airplanes within the SPD group 
given a commanded climb rate of 500 ft/min, a gi irn i t of 1 .25 g’s, and r of 60 seconds. The 
intruder is flying at a vertical rate of either 2750 ft/min or -2750 ft/min. This plot clearly shows the 
impact of intruder velocity, azimuth, and vertical rate on CPA. Comparing the plots of the two 
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Figure 22: SPD Group; r= 60 seconds; gu m i t = 1-25 g’s; Commanded Climb Rate = 500 ft/min; 

Intruder Vertical Rate = ±2750; Segregated by Intruder Azimuth and Intruder Velocity; 
Intruder Velocity Shown by Color. 


intruder vertical rates containing head-on intruders (/ 3 = 0°) reveals the descending intruder 
enables more separation than the climbing intruder. This observation carries over to the 
encounters involving an intruder approaching from p = 45°. For an intruder approaching abeam 
to the UAS (P = 90°), the vertical rate of the intruder appears to have very little effect on CPA. 

While these differences produced “fingers” in the data patterns of the full SPD group 
dataset, the differences in both intruder climb rate and /? are small. Thus, after inspection they 
may be ignored in global analysis of the SPD group, and it is likely that a response surface 
model will be accurate enough to be useful. 

3. 2.3. 2 Response Analysis 

As shown in the previous section, the effects of multiple parameters used to define the 
encounter between the UAS and intruder resulted in small variations in the CPA response. It is 
therefore expected that the overall design space is readily predictable. Figure 23 shows the 
initial r used to define the encounter versus the actual time to CPA. The change in encounter 
time increases as r increases. This follows from the discussion in Section 3. 2. 3. 1.1 which 
showed that the impact of gu mit on CPA increases with t. In order to identify the cause of the 
change in encounter time, Figure 24 shows only the 120-second encounters. The color 
represents the commanded climb rate. From the figure, the larger commanded climb rates result 
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in further deviation from the predicted initial r. Again this follows from discussion in Section 
3.2.3.I.I. 

Figure 25 shows the same design space of 120-second encounters as Figure 24, but 
separates and colors the plot by intruder vertical rate. The figure labels the initial r value on the 
x-axis using a vertical dashed line. From the figure, intruders descending at -2750 and -5000 
ft/min reduce the time to CPA. All other intruder vertical rates tend to increase the time to CPA 
and thus provide additional maneuvering time. 
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Figure 23: SPD Group; t vs Time of CPA. 
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Figure 24: SPD Group; tvs Time of CPA; r = 120 seconds. 
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Tau (sec) vs. Time of CPA (sec) 
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Figure 25: SPD Group; Initial t vs Time of CPA; r = 120 seconds; Segregated by Intruder 
Vertical Rate. 


3. 2.3. 3 Parametric Equation 

The SPD group consists of cases where airplanes have enough energy to maintain the 
commanded climb rate without sacrificing forward speed. As shown in detail in the previous sub- 
sections, these airplanes have predictable climb performance defined by the airplane’s g t imit , 
the z, and the commanded climb rate. Applying response surface methodology using the JMP® 
software, the authors created an equation capable of predicting the achievable CPA for the 
configurations and encounters tested with sufficient accuracy to be useful. The result is 
Equation 32. The equation is of second order containing three primary independent variables 
and their associated interactions. Entering a matrix of the three independent parameters into the 
equations results in the data presented in Figure 26. 

CPA = - 818.76 + 14 * t + 0.93224 * h + 12.4727 * gumit actucd + (h- 837 . 45 ) * 

{ 0.0171 * (t - 57 . 543 ) + 0.03456 * (gumit actual ~ 1 . 8723 )} Equation 32 

Figure 26 shows the relationship between the commanded climb rate and achievable CPA 
for airplanes in the SPD Group. The CPA range is limited to 1500 ft to focus on the area of 
interest. The figure uses color to represent various r values. As r increases, an airplane’s 
achievable separation increases for a constant climb rate. Each r value is represented 
graphically by three lines with various gu m u restrictions: a lower line of 1 .25 g’s, a white dotted 
line of 2 g’s, and an upper line of 3 g’s. The ‘coning’ evident with constant r represents the 
effect of gu m i t on the response. 

The same set of results is presented in a different relationship in Figure 27. This figure 
shows the CPA response as a function of r. The effect of gu m i t is represented as the constant 
thickness within a commanded climb rate and variations in thickness between climb rate values. 
For example, the 500-ft/min commanded climb rate shows that gu mit has very little effect on the 
achievable separation for such a benign maneuver. Table 3 shows that the time required to 
pitch to achieve the commanded maneuver is very small compared to the entire maneuver time 


36 


of the simulation. The fact that the effect of gu m i t is constant for a given maneuver is also 
shown in Figure 26 as the constant height for each r for a single x-axis value (commanded 
climb rate). 
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Figure 26: Application of SPD Group Parametric Equation; CPA vs Commanded Climb Rate; t 
S hown by Color. 
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Figure 27: Application of SPD Group Parametric Equation; CPA vs r; Commanded Climb Rate 
Shown by Color. 


37 


The results shown in this section provide a useful relationship between CPA and three key 
parameters: t, commanded climb rate, and gu m i t . These results are valid for any UAS and 
encounter in the study as long as the UAS/encounter combination is in the SPD group. Thus, a 
UAS that can sustain a commanded climb rate in Figure 26 or Figure 27 from its current flight 
condition can use the predictive equation to determine when it must initiate the climb to achieve 
a particular CPA. The caveat is that results are for a precise and immediate maneuver initiation 
and execution. Also key is that many UAS will not be in the SPD group if the h command is very 
large. Even many single-engine general aviation aircraft cannot sustain climb rates above 500 
ft/min without sacrificing significant speed. This brings us to the next and much more difficult 
performance group. 

3.2.4 Insufficient Power Differential Group 


The airplanes that occupy this performance group attempt to climb at the commanded rate, 
but are forced to sacrifice forward velocity to do so. These airplanes may be required to adjust 
their climb rate to V BA at some time during the encounter. Depending on the specific airplane, 
initial condition, and encounter, they may continue to lose speed throughout the encounter 
without reaching V BA - The performance of individual airplanes within this group is very difficult to 
predict without running the time history. A goal of this study was to model them using multi- 
dimensional modeling techniques based on numerous results from the 2PAIRS simulation and 
to gain insight on relationships between various airplanes in this group and SAA system 
requirements. It was known that this was a very difficult task and no suitable response surface 
model has been developed to date. However, a technique for narrowing down the minimum 
requirements and understanding the details of the parameter relationships within this group has 
been developed and is presented. 

Figure 28 shows all of the CPA responses for the I PD group as a function of specific excess 
power with the result color coded by commanded climb rate. Note that commanded does not 
equal achieved or sustained rate for this group like it did in the previous group. Analysis of this 
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Figure 28: Design Space with only IPD Group; Commanded Climb Rate Shown by Color. 
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performance group focuses on encounters which result in CPA values below 2000 ft. Figure 29 
shows the CPA responses for the I PD group limited to 2000 ft. As P Si increases, the 
performance of the aircraft approaches that of the SPD group, which results in the horizontal 
clusters of points. Figure 30 shows the same design space segregated by t. It is difficult to 
glean conclusive information from either of these plots. The sections that follow dig deeper to 
gain insight and draw some conclusions. 
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Figure 29: Design Space with CPA Range with only IPD Group; Commanded Climb Rate 
Shown by Color. 


CPA (ft) vs. Specific Excess Power (ft/s) 



°S?§ °S?§ °S?§ °S?§ 


Commanded Climb Rate (fpm) 

1 500.0 

1000.0 

1500.0 

2000.0 

2500.0 
30000 
I 3500.0 
H 4000.0 
H- 4500.0 

5000.0 


Specific Excess Power (ft/s) 


Figure 30: Design Space with CPA Range with only IPD Group; Segregated by Time-to-go; 
Commanded Climb Rate Shown by Color. 
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3. 2.4.1 Analysis of Parameter Effects on Climb Maneuver for I PD Performance Group 


The IPD group is characterized by airplanes that are unable to maintain the commanded 
climb rate without sacrificing speed. Within this group a hierarchy of parameter effects on CPA 
exists. The most influential parameters are discussed in this section. Many parameter effects for 
this performance group are very similar to the effects on the SPD group, as described earlier. 

3. 2. 4. 1.1 Effect of Intruder Vertical Rate 


Unlike the SPD group, some intruder elevation angles can cause large variations in the 
achieved CPA. Figure 31 shows the Insufficient Power Group design space segregated by 
intruder vertical rate with r shown by color. Figure 31 reveals an area of high t values with near- 
zero CPA for encounters with an intruder climbing at a high rate. These intruders are climbing 
toward the UAS altitude when the simulation starts and will reach it at the predicted x. Once the 
simulation starts, both aircraft are climbing and the intruder may be climbing faster than the 
UAS. Additionally, in this performance group, the UAS is slowing down. Thus, with the right 
combination of intruder vertical rate, time-to-go, and ownship performance limitations, the two 
airplanes may collide during the maneuver. 
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Figure 31: Design Space with CPA Range with only IPD Group; Time-to-go shown by 
Color; Segregated by Intruder Vertical Rate. 


As an example of a collision caused by this effect, Figure 32 shows the trajectory of the UAS 
and intruder during such an encounter. In the example, the UAS is initially flying at 176 KTAS 
and given a 5000-ft/min commanded climb rate. The intruder is flying at 250 KCAS and climbing 
at 5000 ft/min. The UAS begins the simulation at 6666.7 ft MSL with 80 seconds until collision. If 
either airplane were given a different commanded climb rate, the collision would not happen. 
Obviously, a different avoidance maneuver would be more appropriate for this encounter. Note 
that predicting the collision would likely require a trajectory generator currently. However, use of 
results from this data may be useful in producing guidelines for predicting the bad outcome with 
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less computational expense given the current flight state, performance information for the UAS, 
and the intruder’s trajectory. 



Figure 32: Example GUI Output Showing Example of Collision after Maneuver. 


3. 2. 4. 1.2 Effect of g Urnit 

Figure 31 depicts cluster pairs of maximum achievable CPA values for each x. Figure 33 
reveals the cause of these cluster pairs is common gi irn i t values in the data. Figure 34 is a 
closer look showing only cases with r = 10 seconds with gi irn i t shown by color. It is important to 
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Figure 33: Design Space with CPA Range with only IPD Group; gu mit Shown by Color. 
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note that the resulting analysis of r = 10 seconds extrapolates to the larger r values. Since only 
two values were used to span the gu m i t range of interest, there are only two colors in Figure 34. 
From the figure, there are several airplanes with high values of gi lmit (i.e., red in the figure) that 
seem to follow the trend of lower gi lmlt airplanes (i.e., blue in the figure). 
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Figure 34: Design Space of r = 10 seconds; with CPA Range with only IPD Group; gu m i t 
Shown by Color. 

To address this, Figure 35 shows the same plot with the color legend now showing the 
aerodynamic gu mit the airplanes can achieve rather than the commanded gu m i t , as described 
in Section 2.2.7. This created a third value of gu mit equal to 2.2727 g’s. The total of only three 
values of gu m i t attainable is due to the coarseness of the current design space. Due to the 
multi-dimensionality of the design, most parameters were limited to only 3-5 step changes within 
their respective ranges, as shown in Table 2. Thus, there were only 3 steps of % C, . which 
dictates the initial velocity and flight condition, for all airplane configurations. The values used 
were 8%, 44% and 80%. The lowest value was an assumed value for cruise condition, the 
highest for loiter, and 44% is the midpoint. From Equation 6, these % C Lmnv values correspond 

1 1 L CtA 

to g u m i t values of 12.5, 2.27, and 1.25, respectively. The maximum gu mit commanded in the 
test matrix is 3 g’s, and therefore, the maximum limit is well below the achievable 12.5 g’s for 
some configurations. However, when the 3-g command cannot be achieved the aircraft will be 
aerodynamically limited to either 1 .25 or 2.27 g. 

Even with the separation of attainable gu mit values at a single t, there still exist multiple 
groups of each gumu- This is due to the commanded climb rate and is discussed further in the 
following section. 
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CPA (ft) vs. Specific Excess Power (ft/s) 
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Figure 35: Design Space of r = 10 seconds; with CPA Range with only IPD Group; gu mit 
Actual Shown by Color. 


3. 2. 4. 1.3 Effect of Commanded Climb Rate 

As previously mentioned, there are multiple clusters of points in Figure 35 for a single 
attainable g ii7nit and t combination. These clusters can be shown, as in Figure 36, to be due to 
various commanded climb rate values. With the commanded climb rates separated, there is 
only a single cluster for each attainable gu m i t and r combination. 

Figure 36 clearly shows the quantity of runs for each commanded climb rate. An artifact of 
the performance group, there are fewer airplanes with a 500-ft/min climb rate because more 
airplanes are able to sustain the low climb rate, which places those configurations in the SPD 
group. 
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Figure 36: Design Space of r = 10 seconds; with CPA Range with only IPD Group; g U mit 
Actual Shown by Color; Segregated by Commanded Climb Rate. 

3.2. 4. 1.4 Effect of AV and P s . 

As previously mentioned, a key indicator of climb performance is P Sr Figure 37 shows clear 
minimum values of P Si that will assure a minimum CPA value for each x. Essentially the large 
number of points form a lower right boundary with all larger P Si values yielding a minimum CPA 
value or better. Values of P Si lower than a selected value will yield some CPAs larger and others 
smaller. Thus, a minimum P Si is an indicator that can assure a minimum separation or be used 
to determine a minimum x to initiate a maneuver for a desired separation. For instance, with x = 
30 seconds, assuring a CPA greater than VFR-IFR separation of 500 ft is accomplished with P Si 
values of 28 ft/sec or greater without regard to other parameters. Thus, if the airplane has a 
larger value of P Si , the airplane will be able to achieve the separation of 500 ft for x > 30 
seconds. Some airplanes with lower P Si can also generate CPA > 500 ft, but it is not assured. 
This single parameter is too limiting on the design space because other combinations of 
parameters can also meet the same requirement. One thing the aircraft with lower P Si values 
have in common is speed available (AK) to sacrifice for altitude. The interaction between P Si and 
AV was the apparent area to explore. 

The P Si value necessary to achieve a 500-ft separation is shown to be applicable to 
airplanes with very little velocity to sacrifice, meaning the initial velocity is near the limiting 
velocity, V BA . If an airplane has a high AV, a smaller P Si is required to achieve the same CPA. 
This is shown graphically in Figure 37 as red and grey points above and to the left of the blue 
line. 

As the value of AV increases, the CPA response increases. This relationship indicates that 
the amount of speed an airplane is able to sacrifice directly affects the achievable CPA. Figure 
38 shows a subset of the previous plot, focusing on P Si less than 30 ft/sec, x = 30 seconds and 
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segregates the data by h command . From this figure, there exists a P Si value that can guarantee a 
specified separation for a commanded climb rate. For instance, if an airplane attempts a 2000 
ft/min climb, a P Sl > 26 ft/s will guarantee 500 ft of separation from the intruder. Further, there 
are different P Si values for various intruder azimuths, as shown in Figure 39. This figure shows 
airplanes given a 2000 ft/min commanded climb rate segregated by intruder azimuth. From the 
figure, the ownship requires roughly 26 ft/sec, 18 ft/sec and 18 ft/sec for intruder azimuths of 0°, 
45°, and 90°, respectively. There is a large deviation in required P Si between encounters of 0° 
and 45°, while there is little to no deviation between encounters of 45° to 90°. Thus, future work 
should address additional increments between 0° and 45°. 
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Figure 37: Design Space with CPA Range with only IPD Group; AV Shown by Color; 
Segregated by r; Limited to P s < 30 ft/sec. 
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Figure 38: Design Space with CPA Range with only IPD Group; AV Shown by Color; t = 30 
seconds; Limited to P s < 30 ft/sec; Segregated by h. 
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Figure 39: Design Space with CPA Range with only IPD Group; AV Shown by Color; r = 
30 seconds; Limited to P s < 30 ft/sec; Commanded Climb Rate = 2000 fpm; 
Segregated by Intruder Azimuth. 
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Focusing on the head-on encounter from Figure 39, a P s > 26 ft/sec will guarantee a 
separation of 500 ft; however, there exist airplanes with lower P s . values that are able to achieve 
more than 500 ft of separation. To examine this further, Figure 40 shows the airplanes with 
Pc. less than 26 ft/sec. The figure relates CPA to the amount of forward speed the airplanes 
have to sacrifice for vertical speed. From the figure, if the airplane has P s less than 26 ft/sec, a 
500-ft separation can be achieved if the airplane has more than 55 KTAS to sacrifice, or 55 
KTAS between initial speed and V B a- 
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Figure 40: Design Space with CPA Range with only IPD Group; AV Shown by Color; r = 30 
seconds; Limited to P s < 26 ft/sec; Commanded Climb Rate = 2000 fpm; Flead-on 
Encounter. 


The described performance parameters are capable of providing insight into the CPA a 
specified airplane configuration can achieve during a climb avoidance maneuver without 
sufficient power differential to maintain the commanded rate. Provided enough time, excess 
power or the ability to sacrifice forward speed for vertical, an airplane can predict an achievable 
separation that assures safe separation from an intruder on a collision course. 
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4 Results and Analysis of Descent Maneuver 

This section presents the data from the vertical descent maneuvers along with analysis of the 
results. For reference, Table 4 shows the design space overview for descent maneuvers. 
Descriptions of key influential performance parameters are included as they are presented. 


Table 4: Simulation Parameters Relevant to Descent Maneuver. 


Parameter 

Parameter Symbol 

Values 

Wing loading (lbs/ft 2 ) 

w 

Y 

1,9, 17, 25 

Thrust-to-weight ratio 

T 

w 

0.1, 0.2, 0.3, 0.4 

Aspect ratio 

AR 

5, 16.67, 28.3, 40 

Maximum coefficient of lift 

L max 

0.6, 1.3, 2 

Initial percent of maximum 
coefficient of lift 

°^ C ^max 

8, 44, 80 

Maximum lift-to-drag ratio 

L 

Dmax 

5, 20, 35, 50 

Load factor Limit 

dlimit 

1.25, 3 

Altitude (ft MSL) 

h 

1000, 6666.7, 12333, 18000 

Commanded climb rate 
(ft/min) 

h command 

500, 2000, 3500, 5000 

Intruder velocity (KCAS) 

Vintr 

50, 150, 250 

Intruder vertical rate 

hintr 

-5000, -2750, -500, 0, 500, 2750, 
5000 

Intruder azimuth (deg) 

P 

0, 45, 90 

Time-to-go (seconds) 

T 

10, 20, 30, ... 120 

Velocity Not to Exceed 

K NE 

1.15* V cruise 


4.1 Descent Maneuver Performance Parameters 

The descent maneuver shares some performance parameters with the climb maneuver 
though they are generally associated with opposite signs. It has some additional considerations 
as well. The most significant difference is obviously that the airplane tends to speed up when 
descending vs. tending to lose speed in a climb. This section describes key performance 
parameters that limit descent maneuvers. 

4.1 .1 Specific Excess Power (P s ) 

The primary performance parameter for the climb maneuver was Ps imax , as defined in 
Equation 25. This equation defines the airplane’s specific excess power as a function of the 
maximum thrust available. The power parameter of interest for the vertical descent maneuver is 
associated with minimum thrust: 


Pc. . 

min 



Equation 33 


For this study, the minimum thrust of all airplane configurations is zero. Consequently, 
Equation 33 becomes: 
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Pc. . =(- — }*Vj. Equation 34 

In the climb maneuver this parameter determined how fast the aircraft could climb without losing 
speed. Here the parameter determines how fast the aircraft can descend without gaining speed 
as well as how rapidly and how much speed it will gain if constant speed is not possible. Akin to 
the climb maneuver, there is a maximum velocity an aircraft cannot exceed, V NE , which is 
discussed in more detail below. For the discussion of vertical descent performance, all 
references to Pc. refer to the Pc. . as calculated in Equation 34. 

4.1.2 Potential Change in Velocity (AV) 

An aircraft can have a number of maximum speeds; maximum turbulence penetration 
speed, maximum flap extend speed, and velocity to never exceed (V NE ) are key ones. V NE is the 
only one of interest to this effort. It is related primarily to structural limitations normally, but when 
it is a maximum Mach number it may be related to stability and control. Either way, it can limit 
the ability to descend at high rates or at high rates for extended periods of time. Like the 
potential change in velocity for the climb, the descent maneuver has the same parameter but it 
is related to Vne- The amount of speed the airplane can gain affects the achievable separation. 
This parameter is a measure of the speed a given airplane can gain before reaching the Vne- 
The equation for this parameter is: 

AV = V NE - Vi. Equation 35 

This parameter in association with Ps imin is directly related to the rate and duration of 
descent an airplane can maintain to generate separation from an intruder. 

4.1.3 Difference between Commanded Descent Rate and Accelerating Descent Rate 

(AflAccelerate) 


For a given airplane configuration there is a maximum descent angle and rate above which 
the airplane cannot avoid accelerating from its initial speed even at idle power. If the airplane is 
commanded to descend at a smaller angle, the airplane will be able to descend without 
increasing velocity. The difference between the commanded descent rate and this sustainable 
rate will provide insight into the rate of acceleration the airplane configuration will experience 
during the maneuver. This performance parameter can be expressed mathematically as: 

Ah ac celerate ~ h command ~ h accelerate ~ h command ~ Equation 36 

This parameter is similar to, and has a similar impact as the Ah Sustainab i e parameter 
discussed in the vertical climb maneuver section. 

4.2 Performance Grouping of Descent Maneuver Data 

Initial analysis revealed multiple performance groups independent of the airplane design 
parameters. Several different ways of grouping the descent maneuver data were considered. 
Initial response surface models were fit to preliminary data and residual analysis was done. 
These groups were revealed to be negatively affecting the accuracy of the model fit. Using 
JMP® data exploration tools and 2PAIRS GUI functionality, individual airplane configurations 
within each group were analyzed. 

Three predictable performance groups presenting three different characteristic flight 
performances during the commanded maneuver were created. The process for discovering and 
segregating the three performance groups was similar to the process outlined in Section 3.2, 
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and the resulting groups are analogous to three of the groups identified for the climb maneuver: 
a Sufficient Power Differential (SPD) group, an Insufficient Power Differential (IPD) group, and 
an Insufficient Time to Pitch (ITP) group. The three performance groups are detailed in the 
following subsections. 

Like the climb maneuvers, an airplane configuration alone does not set the performance 
group. The combination of airplane design parameters, commanded maneuver rate, initial flight 
condition, and encounter parameters are all required to determine performance group 
association as it relates to the SAA system requirements. Additionally, small deviations within 
an airplane’s design parameters or the commanded maneuver may result in that configuration 
changing performance group. Association with each of the performance groups is predictable 
prior to the initiation of the 2PAIRS simulation runs using the aircraft design, encounter, and 
maneuver parameters. 

4.2.1 Insufficient Time to Pitch Group 

Similar to the ITP climb maneuver performance group, the airplane configurations that fall 
within this grouping are unable to reach the commanded descent rate during the given 
encounters. This grouping is dominated by fast-moving airplanes with low gi imit involved in 
encounters with short t, generally 10 seconds in the data evaluated to date. These same 
airplane configurations are found in the other performance groups with higher t values. 

4.2.2 Insufficient Power Differential Group 

For the cases in this group, the airplane configurations attempt to fly at the given 
commanded descent rate and gain speed during the maneuver because the commanded 
descent rate is greater than their h accelerate . Due to the acceleration, these airplanes may be 
required to adjust their descent rate to a sustainable rate if V NE is reached. Given a separation 
threshold of interest, an analysis approach analogous to the IPD group done in the climb 
avoidance maneuver analysis can be used for airplanes within this group. Thus, once a 
specified separation threshold is established, analysis can be completed similar to that outlined 
in Section 3.2.4. 

Within this performance group, there exist several characteristic maneuver scenarios. Given 
adequate time, the airplanes will either: 1) reach terminal velocity for the commanded descent 
rate and maintain the commanded descent rate throughout the simulation, or 2) reach their 
specified V NE and be forced to adjust to a sustainable descent rate. Which of these maneuver 
scenarios occurs is highly dependent on the airplane’s drag configuration and the commanded 
descent rate. Without adequate time, the airplane will continue to accelerate throughout the 
simulation run reaching CPA before it reaches V NE . Currently, it is difficult to predict which of 
these scenarios will occur as it is highly dependent on the airplane configuration, commanded 
maneuver, and t value of the encounter. 

4.2.3 Sufficient Power Differential Group 

For the cases that occupy this performance group, the airplane configurations do not 
accelerate during the commanded descent maneuver because they are given a commanded 
descent rate that is less than their h acce i erate . Therefore, the Ah accelerate value for cases within 
this performance group is positive. This group is directly analogous to the climb SPD group 
described in Section 3.2.3. 

The commanded descent rate is very influential on when airplane configurations fall within 
this group. The higher the commanded descent rate, the greater the influence of the gravitation 
force on the acceleration of the aircraft. Thus, there exist more airplane configurations with a 


50 



500-ft/min commanded descent rate than those with a 5000-ft/min descent rate within this 
performance group. 

The lack of acceleration results in highly predictable separation performance for these 
airplane configurations. Using JMP® software, a parametric equation was developed to predict 
CPA as a function of x, commanded descent rate, and gu mit . The equation is: 

CPA = -1170.71 + 17.62 * r + 1.06 * h + 19.08 * g limit + (h- 1058.3) * [. 0167 * (t - 65) + 
.0365 * ( gu m i t - 1.96)]. Equation 37 


Figure 41 and Figure 42 show the application of the parametric equation to a matrix of the 
three parameters. The use of these figures enables quick prediction of achievable separation. 
Figure 41 shows CPA versus the commanded descent rate, while Figure 42 shows CPA versus 
t. The gu m i t effect on these figures is similar to that described for the climb SPD Group 
described in Section 3. 2. 3. 3 

From Figure 41 , if an aircraft is able to detect the intruder with x = 20 seconds, with the 
knowledge that it can perform a 3000-ft/min descent without accelerating, the airplane can 
achieve between roughly 800 and 1000 ft of separation depending on the gu m i t . If further 
separation is required, the airplane must either have a longer detect time to CPA or have an 
increased sustainable descent rate. A similar application can be used for Figure 42. 
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Figure 41: Application of SPD Group Parametric Equation; CPA vs Commanded 
Descent Rate; r Shown by Color. 
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Figure 42: Application of SPD Group Parametric Equation; CPA vs r; 
Commanded Descent Rate Shown by Color. 


52 


5 Results and Analysis of Level Turn Maneuver 


This section presents the data from simulation runs of the level turn maneuver along with 
analysis of the results including exploration of various parameter effects on CPA. For reference 
during the following discussions, the design and encounter parameters relevant to the descent 
maneuver simulations are listed in Table 5 showing the discrete values simulated. These 
parameters were described above in Sections 0 and 0. 


Table 5: Simulation Parameters Relevant to Level Turn Maneuver 


Parameter 

Parameter Symbol 

Values 

Wing loading (lbs/ft 2 ) 

w 

~s 

1, 9, 17, 25 

Thrust-to-weight ratio 

T 

W 

0.1, 0.2, 0.3, 0.4 

Aspect ratio 

AR 

5, 16.67, 28.3,40 

Maximum coefficient of lift 

^ L max 

0.6, 1.3, 2 

Initial percent of maximum 
coefficient of lift 

%C Lmax 

8, 44, 80 

Maximum lift-to-drag ratio 

L 

D max 

5, 25, 35, 50 

Load factor Limit 

9 limit 

1.5 

Altitude (ft MSL) 

h 

1000, 6666.7, 12333, 18000 

Commanded climb rate (ft/min) 

h command 

500, 2000, 3500, 5000 

Intruder velocity (KCAS) 

Vintr 

50, 150, 250 

Intruder vertical rate 

hintr 

-5000, -500, 0, 500, 5000 

Intruder azimuth (deg) 

P 

-90, -45, 0, 45, 90 

Time-to-go (seconds) 

T 

10, 20, 30, ... 120 

Turn Rate Command (deg/sec) 

U* command 

1,3, 5, 10 

Roll Rate (deg/sec) 

0 

30 

Heading Change (deg) 

L 10 

90 


Attempts to develop a response surface model have not resulted in adequately accurate 
representations to date. Results from the large number of runs made are presented and 
discussed. Useful conclusions are drawn for the segment of airplanes evaluated. Analysis only 
includes airplanes doing maneuvers that are not limited by their g t imit and are able to perform 
the commanded maneuver without changing speed. This approach eliminates artificial limits that 
impact evaluation of trends. The gumit enforced for the results presented is 1 .5 g, which is 
equivalent to a 48° bank-angle limit. The encounter scenarios and aircraft configurations are not 
exhaustive, but are representative of the key behaviors of interest. 

The level turn avoidance maneuver is the only maneuver studied that is not used by TCAS 
to prevent collisions. It is commonly used by pilots to remain well clear of other aircraft. Also, 
small deviations are generally less disruptive to the NAS. Thus, although it is less likely to be 
the better choice for collision avoidance than a vertical maneuver, it is usually the first choice for 
remaining well clear. 

As detailed in Section 2.5, the acceptable horizontal separation is difficult to discern from 
current separation standards. The section concludes that the range of interest for horizontal 
CPA results encompasses 500 ft and 3 nm of horizontal separation. 
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5.1 Parameter Effects on Level Turn Maneuver 


There are several possible limitations to the turn performance of an airplane. For a given 
airplane configuration, the turn rate can be expressed as: 


ip = ^* Vnf - 1 = g * 


V 


-*p*C L 

W. 

s 



Equation 38 


As the equation shows, for a fixed — , and p, there is a fixed V and C L combination that 

must be met for the flight condition to be valid. 


Similarly, the turn radius can be expressed as: 



p*g*C L 



V 


Equation 39 


From Equation 38, the airplane configuration with fixed — and large Q, and low 

altitude (high p) lead to higher turn rates. The encounter scenario, which includes the intruder 
parameters, altitude, and t, also impacts the separation performance for a given scenario. 
Examples of these effects are shown in the sub-sections below. 


5.1 .1 Effect of Altitude 


Equation 39 shows the relationship between the turn rate, turn radius, and velocity. As 
altitude increases, the true airspeed (KTAS) of the UAS increases for a constant calibrated 
airspeed (KCAS). Thus, for a constant turn rate, as the altitude increases, the radius of the turn 
increases. Figure 43 shows this effect on the CPA. The ownship and intruder have the same 
KCAS for both results. The UAS is given a constant commanded turn rate for both encounters, 
but the altitude varies from 1000 ft MSL (blue line) to 18,000 ft MSL (red line). The lower-altitude 
UAS makes a tighter turn and initially generates a larger CPA. However, the faster aircraft may 
ultimately generate a larger CPA over longer t. Note that aerodynamic performance limitations 



Figure 43: Level Turn Maneuver; Altitude Effect on Response. 
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that drive turn rate are associated with KCAS while this demonstrates that the resulting turn 
performance (radius) is linked to KTAS. 

5.1 .2 Effect of Intruder Azimuth 

The intruder’s azimuth intuitively was expected to have a large effect on the achievable 
separation. For the level turn maneuvers in this study, the intruder may approach perpendicular 
to the ownship’s trajectory from either side. With the ownship only turning right, if the intruder is 
initially right of the ownship’s trajectory, the UAS may turn directly into the intruder’s constant 
flight path. Trivially, this may create a smaller CPA value than if the ownship were to turn 
directly away from the intruder. The only encounter scenario where this may be a recommended 
avoidance maneuver is if the ownship can turn and fly parallel to the intruder with a separation 
greater than the required threshold. 

Conversely, if the ownship turns away from an intruder that is approaching at a 90° azimuth, 
the UAS may gain separation from the intruder depending on the relationship between the 
ownship’s speed and the intruder’s speed. For example: if the UAS is flying faster than the 
intruder, it will turn to the commanded 90° heading change and will create further separation 
from the intruder. 

Figure 44 shows a single UAS configuration involved in an encounter with intruders 
approaching from an array of azimuths. Each intruder is slower than the UAS and all have the 
same speed. Thus, intruder azimuth is the only parameter that varies between the cases shown 
in the figure. The intruders’ trajectories are represented by the colored dashed lines with the 
initial positions marked by a star. The CPA for both the UAS and intruder is represented as a 
colored dot. If no avoidance maneuver were performed, a direct collision (CPA = 0 ft) would 
occur where the intruder trajectories intersect. 



Figure 44: Level Turn Maneuver; Intruder Azimuth Effect. 

The UAS is able to achieve a separation of greater than 0.5 nm regardless of the angle from 
which the intruder approaches, for the given encounter parameters. Although the CPA occurs 
before the initial r value, a greater amount of separation occurs and a direct collision is avoided. 
This shows that intruder azimuth does not have a large effect when the intruder is significantly 


55 


slower than the UAS. The next section discusses other effects of the speed ratio between the 
UAS and intruder. 


5.1 .3 Effect of Intruder Speed 


The magnitude of the intruder’s speed effectively changes the closure rate between the 
ownship and intruder for a given UAS configuration. This effect is especially influential if the 
ownship is turning away from the intruder. Figure 45 shows the same encounter as Figure 44, 
but the intruder’s speed is increased to 198 KTAS, significantly faster than the UAS. Comparing 
Figure 44 and Figure 45, one significant difference is the time of CPA for the intruders 
approaching from the right (-90 and -45°). For the -90° intruder, once the UAS reaches a 90° 
heading change, the two aircraft are flying parallel to each other. If CPA has not yet occurred, 
the CPA will occur when the intruder catches the UAS. This is a result of the systematic way the 
maneuvers are being flown with a heading change of 90°. A maneuver algorithm or UAS pilot 
would make a more effective avoidance maneuver based on the applicable mission profile. 
Intruder azimuth had a less consistent effect on the difference between predicted and actual x 
with the faster intruder than the slower intruder. The increase in intruder speed resulted in actual 
x being larger than predicted x for the -90° and -45° azimuth intruders. The remaining azimuths 
resulted in actual x being smaller than predicted x but with very little difference for the head-on 
(0°) case. Thus, turning away from the intruder gives the UAS and its operator more time to 
consider options if the intruder is faster than the UAS. However, depending on x and the turn 
radius, the CPA could become very small. 


XlO 

1 . 5 - 


o 0.5 

c 

cc 

V) 

>- 0 


-0.5- 




V(KTAS)= 142.8434 ,Tau(sec)=30 .Commanded Turn Rate (deg/s)=3 
Intruder: V(knots)= 198 Azimuth (deg)= -90 

T * 


• CPA (fl)=2468 Time of CPA (s)*75 



— Intruder: V( knots)- 198 Azimuth (deg)= -45 



• CPA {ft)=4346 Time of CPA [s)=40 



— Intruder: V(knots)= 198 Azimuth (deg)= 0 



• CPA {ft)=4911 Time of CPA (s)=29 

1 


— Intruder: V(knots)= 198 Azimuth (deg)= 45 
* CPA (ft)- 4242 Time of CPA (s)=26 

Intruder: V( knots)- 198 Azimuth (deg)= 90 
CPA (ft)-2799 Time of CPA [s)*23 


0.5 


15 

X-dtstance (ft) 


2,5 


xlQ 


Figure 45: Level Turn Maneuver; Intruder’s Speed Effect on Various Azimuths. 


5.1.4 Effect of Turn Rate 

Figure 46 throughFigure 49 are contour plots showing achievable CPA for several velocity 
and r combinations. The contour lines in each depict CPAs of 500 ft, 1 nm, 2 nm and 3 nm. For 
this study, a gii mit of 1 .5 was used in the turning maneuver for all airplane configurations. The 
airplanes represented in these contour plots did not reach this limit so they were following the 
commanded turn rate. From Equation 38, the gu m it would be a factor at or above 234 KTAS for 


56 


a 5-deg/sec turn and 120 KTAS for 10-deg/sec turn. The contour plots are not valid for speeds 
above these limiting values and are, thus, shaded out. 


1 deg/sec turn 



Figure 46: Level Turn Maneuver; Head-On Encounter; gu mit < 1.5; h = 18,000 ft MSL; 
V int r = 50 KTAS; Contour Showing Achievable CPA for 1 -deg/sec Turn. 


3 deg/sec turn 



Figure 47: Level Turn Maneuver: Head-On Encounter; gu mit < 1.5; h = 18,000 ft MSL; V intr 
= 50 KTAS; Contour Showing Achievable CPA for 3-deg/sec Turn. 


57 


5 deg/sec turn 



250 


300 


Figure 48: Level Turn Maneuver; Head-On Encounter; g Umlt < 1.5; h = 18,000 ft MSL; V intr 
= 50 KTAS; Contour Showing Achievable CPA for 5-deg/sec Turn. 


1 0 deg/sec turn 



Figure 49: Level Turn Maneuver; Head-On Encounter; gu mit < 1.5; h = 18,000 ft MSL; V intr 
= 50 KTAS; Contour Showing Achievable CPA for 10-deg/sec Turn. 

As an example of how these plots can be used, the dark blue lines added to each graph 
point out the x required for a UAS flying at 100 KTAS to achieve a CPA of 1 nm. Each figure 
provides results from a different commanded turn rate for a head-on (/ 3 = 0°) encounter with a 
50-knot level intruder flying co-altitude at 18,000 ft MSL. Table 6 shows the required r for each 
commanded turn rate and the results are plotted in Figure 50. As expected, as the turn rate 
increases, less time is required to achieve the 1-nm separation, although the relationship 
between turn rate and time required is non-linear. Due to the variation in r required between a 
1 -deg/sec and 5-deg/sec turn rate, further data should be analyzed between these values to 
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determine the knee in the curve. This figure Table 6: Time Required to Achieve Horizontal 
may be used to determine a required turn rate Separation of 1 nm for 100 KTAS 

ability for a set t, or vice versa. For instance, UAS at Various Commanded Turn 

given a t of 40 seconds, the 100 KTAS airplane 
must be able to sustain a turn rate greater than 
3-deg/sec for a 1-nm CPA. 

These contour plots are valid for all airplane 
configurations that sustain the turn through a 90° 
heading change without losing speed and do not 
reach the 1 .5 gu m i t flying at 18,000 ft MSL with a 
50-knot intruder approaching head-on. Varying these constant parameters would result in a 
series of contours applicable to many encounter scenarios. The series of contours in Figures 46 
through 49 can be used to make Figure 50 into a 3-D plot of turn rate, CPA, and t. Additionally, 
contours for other conditions can expand the applicable results. 
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Figure 50: Level Turn Maneuver; Head-On Encounter; gi lmit < 1.5; h = 18,000 ft MSL; V = 
100 KTAS; V intr = 50 KTAS; Required t vs Turn Rate to Achieve 1 nm of 
Separation. 
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6 Concluding Remarks 

This study of the trade space between UAS maneuver performance and SAA system 
performance provides informative results, analysis, and methodologies to help regulators and 
UAS designers develop regulations and capabilities for safe and efficient integration of UAS into 
the NAS. Simulators for assessing the trade space across a broad range of UAS maneuver 
performance capability and encounters were developed for climb, descent, and level turn 
maneuvers. The UAS maneuver performance modeling approach used in developing the 
simulators was described. The simulator has both batch and GUI interfaces supporting: 1) 
evaluations of macro results from large numbers of UAS performance models and encounter 
geometries; and 2) more detailed time history analysis of model/encounter combinations, 
respectively. 

Test matrix parameter ranges for UAS aerodynamic and propulsive parameters, encounter 
model parameters, and intruder states were established based on existing airplanes and near- 
term NAS integration plans as described by RTCA SC-228’s Terms of Reference and Work 
Plan [4,5], The test matrix and reasoning behind it were discussed. 

Discussion of expected closest-point-of-approach ranges of interest was presented and 
used in analyzing results. Results and methodology for analyzing the large volume of data was 
presented for the three maneuvers. Early evaluation made it clear that successful analysis 
would require compartmentalizing UAS performance capabilities into groups. The performance 
capabilities are grouped by UAS aerodynamic/propulsive performance, initial flight condition, 
commanded maneuver rate, and the encounter time to go. An airplane performance model 
alone cannot determine the performance group in which the airplane belongs, because it is also 
dependent on the encounter scenario and commanded maneuver. The Insufficient Power 
Differential (I PD) group, which consists of cases where the UAS has insufficient performance to 
sustain the commanded climb, descent, or turn rate at the initial speed were considered to be of 
most interest for this work. They are also significantly more difficult to analyze as a group than 
the Sufficient Performance Group, where the UAS is able to sustain the commanded maneuver. 
The latter amounts to a geometry problem of lines and curves. Two groups were not analyzed in 
detail because they were deemed unlikely to occur, especially in flight conditions for transit from 
an airport area to Class A or Special Use Airspace. The transit condition is the focus of SC-228 
[4]. 

A limited evaluation was shown using approximate aerodynamic/propulsive performance of 
an existing manned aircraft. This data served as verification and validation of the simulation tool 
and to provide insight on the impact of initial conditions and aerodynamic/propulsive parameters 
on CPA. 

A detailed discussion of the impact of UAS performance on achievable CPA was presented 
for the vertical climb maneuver and to a lesser extent the vertical descent and level turn 
maneuvers. Results demonstrated that specific excess power, P s , of a given aircraft is a key 
parameter in predicting achievable CPA for the vertical maneuvers. Results also show that while 
a minimum P s value can assure a minimum CPA for a given time-to-go, the same or greater 
CPA may be achievable with a lower P s if the aircraft has sufficient speed to trade for altitude 
(kinetic energy to trade for potential energy). The level turn maneuver is much less likely to 
encounter specific excess power limitations for the load factors allowable for most aircraft. Thus, 
level turn analysis focused on aircraft speed and turn rates limited to a 1 .5-g load factor, which 
equates to a maximum of 48° of bank. Future work will attempt to better quantify the 
relationships for the overall test matrix and is expected to be conducted in conjunction with SC- 
228 Work Group 1 - Detect and Avoid. 
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7 Future Work 


The simulators, results, and methodologies presented in this paper attempted to establish a 
basis and framework for determining a minimum UAS performance that can be associated with 
performance of a given sense and avoid system. This was accomplished for climb, descent, and 
level turn maneuvers. However, some of the specific analysis necessarily made assumptions 
about CPA ranges of interest as well as the overall design space. Going forward the authors will 
work with RTCA SC-228 to target analysis and results to the group’s needs. In particular, SC- 
228’s Work Group 1, Detect and Avoid, has sub-groups established to develop system 
requirements. The authors are working with the UAS Performance group and will continue to 
use the tools and methods presented here in the work. 

Additionally, the simplified maneuvers and perfect sensors used in this study yield idealized 
results. Future studies will expand on the lessons learned from this work to study their more 
real-world implications associated with sensor noise and maneuver algorithms. 
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